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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC)  at  the  request  of  AEDC/DOT.  The 
Air  Force  Project  Manager  was  Mr.  D.  A.  Duesterhaus.  The  results  were  obtained  by 
Sverdrup  Technology,  Inc.,  AEDC  Group,  operating  contractor  for  aeropropulsion  testing 
at  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  Project  Number 
D188EW.  The  experimental  analysis  was  performed  from  August  1981  to  October  1983,  and 
the  analytical  analysis  was  completed  in  March  1984.  The  manuscript  was  submitted  for 
publication  on  March  21,  1985. 

Volume  I  contains  the  analysis  and  evaluation  of  the  icing  scaling  equations,  whereas 
Volume  II  is  a  description  and  listing  of  the  computer  code  entitled  “SIMICE.” 
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1.0  INTRODUCTION 


1.1  BACKGROUND 

The  formation  of  ice  on  aircraft  surfaces  occurs  during  flight  through  supercooled 
droplets.  Ice  accretions  on  these  surfaces  usually  degrade  both  aircraft  performance  and 
operational  safety.  For  this  reason,  it  has  become  important  in  the  design  and  certification 
phases  of  system  development  to  evaluate  system  performance  degradation  because  of  icing. 

The  most  acceptable  method  of  evaluating  the  performance  characteristics  of  aircraft 
and  aircraft  components  for  system  certification  would  be  to  conduct  flight  tests  in  natural 
icing  conditions.  Adequate  system  evaluation  through  icing  flight  tests  requires  specific 
weather  conditions,  posing  severe  operational  limitations  because  of  the  low  frequency  of 
their  occurrence.  This  makes  certification  through  flight  testing  time  consuming  and 
expensive  because  large  amounts  of  flying  time  are  required  to  adequately  evaluate  a  system 
over  a  wide  range  of  conditions.  Flight  testing  in  a  simulated  cloud  produced  by  an  icing 
tanker  would  appear  to  be  the  next  most  desirable  method  of  conducting  icing  tests.  In 
practice,  this  method  of  flight  testing  is  severely  limited  by  tanker  operational  limitations, 
lack  of  control  of  atmospheric  parameters  such  as  temperature,  pressure  and  humidity,  and 
the  logistics  of  maintaining  steady  flight  in  a  cloud  produced  by  a  leading  aircraft. 
Therefore,  ground  testing  in  altitude  facilities  has  become  an  accepted  approach  for 
evaluating  aircraft  system  performance  in  icing  conditions. 

The  size  of  many  of  the  aircraft  components  to  be  tested  greatly  limits  the  number  of  test 
facilities  capable  of  conducting  icing  tests.  It  would  be  beneficial  to  the  icing  community  if 
the  feasibility  of  using  scale  models  in  these  tests  could  be  proven.  The  problem  is 
complicated  because  the  icing  process  is  governed  by  the  impingement  of  water  droplets  and 
the  thermodynamics  of  the  freezing  process  on  the  surface.  This  document  gives  the  results 
of  a  study  conducted  at  Arnold  Engineering  Development  Center  (AEDC)  Engine  Test 
Facility  (ETF)  to  experimentally  verify  a  set  of  icing  scaling  equations. 

1.2  SCOPE  OF  INVESTIGATION 

A  research  program  was  conducted  at  the  AEDC/ETF  to  evaluate  the  icing  scaling 
requirements.  The  objectives  of  the  study  were  (1)  to  evaluate  the  equations  governing  the 
icing  process  to  identify  proposed  scaling  parameters,  (2)  to  develop  a  computer  code  to 
solve  the  various  forms  of  the  icing  scaling  equations,  (3)  to  conduct  tests  to  determine 
which,  if  any,  of  the  proposed  methods  produced  scale  ice  accretions,  and  (4)  if  an  accurate 
set  of  icing  scaling  equations  were  found,  to  write  a  final  computer  code  that  could  be  used 
in  icing  tests. 
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The  scaling  verification  tests  were  conducted  using  full-  and  half-scale  circular  cylinders 
and  full-,  1/3- ,  and  1/6-scale  airfoil  sections.  By  applying  a  postulated  scaling  method,  test 
conditions  that  should  produce  scale  ice  accretions  were  calculated  using  the  computer  code 
developed  for  this  study.  The  size  and  shape  of  the  ice  accretions  resulting  from  these  test 
conditions  were  compared  with  full-scale  results  to  determine  the  accuracy  of  the  scaling 
method.  The  objectives  of  the  study  were  met  by  identifying  a  scaling  method  that  produces 
scaled  ice  accretions  over  a  wide  range  of  test  conditions  and  that  can  be  applied  to  a  variety 
of  icing  testing  situations. 

13  OUTLINE  OF  REPORT 

This  report  begins  with  a  definition  of  the  requirements  of  icing  scaling  equations.  These 
requirements  are  then  examined  in  detail  through  the  development  of  mathematical  models 
to  identify  proposed  scaling  parameters.  A  review  of  past  icing  scaling  methods  and  their 
respective  verifications  will  also  be  discussed. 

Following  the  scaling  parameter  development,  the  experimental  verification  procedures 
and  results  will  be  presented.  Application  of  the  scaling  procedure  will  be  illustrated  through 
the  presentation  of  several  scaling  examples.  Certain  limitations  identified  from  the 
experimental  program  will  also  be  discussed. 

2.0  DISCUSSION  OF  ICING  SCALING 


2.1  DEFINITION 

Icing  scaling  is  a  test  procedure  used  to  form  similar  ice  accretions  on  two  geometrically 
similar  objects  under  different  atmospheric  and  meteorological  conditions.  “Similar”  in 
this  case  implies  not  only  identical  geometrical  shapes  but  identical  surface  characteristics, 
such  as  the  type  of  ice  and  roughness.  Icing  is  an  accretion  process  meaning  that  each  layer 
of  ice  is  formed  on  the  previous  layer.  This  implies  that  if  the  ice-accretion  process  can  be 
started  on  a  clean,  model  airfoil  with  surface  characteristics  similar  to  a  full-scale  airfoil,  the 
ice-accretion  process  will  continue  similarly  and  produce  scaled  ice  accretions.  This 
statement  forms  the  major  premise  of  icing  scaling  and  will  be  applied  throughout  the 
development  of  the  scaling  parameters.  To  apply  an  icing  scaling  procedure,  a  set  of 
equations  that  accurately  model  the  beginning  of  the  icing  process  must  be  defined,  and 
parameters  that  relate  the  atmospheric  conditions  of  velocity,  pressure,  and  temperature  and 
the  meteorological  conditions  of  liquid-water  content  (LWC),  droplet  size,  and  icing  time 
must  be  identified. 
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2.2  TYPES  OF  ICING  SCALING 

The  first  and  most  obvious  type  of  icing  scaling  is  that  in  which  the  ice  accretions  are 
formed  on  a  subscale  model  similar  to  those  that  would  be  obtained  at  full-scale.  The  ability 
to  test  subscale  geometries  would  allow  smaller  test  facilities  to  conduct  a  wider  range  of 
icing  tests.  Since  smaller  test  facilities  are  usually  less  expensive  to  operate  than  larger 
facilities,  the  cost  would  be  substantially  reduced. 

A  second  type  of  icing  scaling,  test  parameter  scaling,  is  aimed  at  increasing  the 
simulation  capabilities  of  test  facilities  by  extending  the  Mach  number,  altitude,  and/or 
temperature  ranges.  Icing  scaling  equations  relate  the  test  cell  atmospheric  and 
meteorological  conditions  so  that  similar  ice  accretions  can  be  formed.  These  equations 
indicate  how  ground  test  conditions  could  be  modified  to  avoid  the  operational  limits  of  a 
test  facility,  yet  achieve  the  appropriate  icing  condition.  For  example,  a  facility  may  be 
limited  to  ambient  pressure  but  have  the  need  to  produce  an  ice  accretion  that  would  be 
formed  at  an  altitude  of  10,000  ft.  The  scaling  equations  would  be  applied  to  determine  the 
test  conditions  at  atmospheric  pressure  required  to  simulate  the  icing  that  occurs  at  the 
altitude  flight  condition.  This  type  of  scaling  will  allow  test  facilities  to  offer  a  wider  range  of 
possible  conditions  to  an  aircraft  manufacturer  to  evaluate  their  flight  system. 

The  applications  of  both  types  of  scaling  would  make  them  very  attractive  to  both  the 
test  facility  operator  and  the  aircraft  manufacturer. 

2.3  ICING  SCALING  REQUIREMENTS 

The  requirements  of  icing  scaling  and  the  individual  processes  that  must  be  modeled  can 
be  identified  by  examining  the  ice-accretion  process.  Ice  can  form  on  an  object  during  flight 
through  clouds  containing  supercooled  droplets.  These  droplets  will  follow  trajectories 
determined  by  the  forces  acting  upon  them  and  will  either  strike  or  miss  the  object.  The 
impacting  droplets  will  either  freeze  in  place,  run  back  along  the  surface  and  freeze  in 
another  location,  or  run  back  and  be  blown  off  the  surface.  The  shape  of  the  resulting  ice 
accretion  will  be  determined  by  the  distribution  of  the  mass  of  water  striking  the  airfoil,  the 
locations  at  which  that  water  freezes,  and  the  type  of  ice  formed.  From  this  elementary 
description  of  aircraft  icing,  the  parts  of  the  ice-accretion  process  that  must  be  considered  in 
an  icing  scaling  analysis  can  be  identified.  They  are  (1)  the  flow  field  about  the  body,  (2)  the 
droplet  trajectory  and  impingement  characteristics,  (3)  the  total  mass  of  water  impinging  on 
the  surface,  and  (4)  the  thermodynamics  of  the  freezing  process.  Each  of  these  areas  will  be 
evaluated  in  subsequent  sections  to  identify  the  scaling  parameters  used  in  this  study. 
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3.0  IDENTIFICATION  OF  THE  SCALING  PARAMETERS 
3.1  FLOW  FIELD  SCALING  REQUIREMENTS 

First,  for  the  geometrically  similar  surface  conditions  required  to  produce  similar  ice 
accretions,  it  is  assumed  that  the  flow  field  over  the  bodies,  i.e.  normalized  pressure  and 
temperature  distributions,  must  be  similar.  This  can  only  be  assured  if  the  two  bodies  have 
scaled  geometries,  and  the  viscous  effects  are  geometrically  scaled.  The  exact  duplication 
of  these  characteristics,  obtained  by  application  of  the  aerodynamic  similitude  laws,  places 
stringent  restrictions  on  the  test  conditions.  In  general,  these  requirements  are  not  compatible 
with  the  other  icing  scaling  parameters  to  be  discussed  in  the  following  sections.  An  additional 
assumption  must,  therefore,  be  made  to  relax  these  requirements. 

On  an  unheated  body,  an  ice  accretion  will  generally  form  only  in  the  region  of  droplet 
impingment.  It  can  therefore  be  assumed  that  the  normalized  pressure  and  temperature 
distributions  need  to  be  similar  only  in  this  region.  Except  for  extreme  angles  of  attack  or 
complex  geometries,  the  droplets  generally  strike  the  body  in  the  vicinity  of  the  stagnation 
point  where  the  boundary  layer  is  relatively  thin.  While  there  are  Reynolds  number  effects 
on  the  boundary  layer  even  in  the  stagnation  region,  the  effect  on  the  ice  accretion  could  not 
.b5...quantified  a.Prif>r^’  and  therefore,  it  was  assumed  that  viscous  effects  are  adequately 
scaled  simply  by  requiring  that  the  bodies  are  geometrically  scaled.  The  applicability  of  this 
assumption  was  not  explicitly  evaluated  at  the  beginning  of  this  study  because  its  validity 
was  determined  by  comparing  the  full-  and  subscale  ice  accretions.  In  general,  the  applicability 
of  this  assumption  for  complex  test  articles  and  test  conditions  which  may  include  high  angles 
of  attack  should  be  proven  either  experimentally  or  analytically  prior  to  conducting  the  icing 
test. 


The  need  for  similar  flow  fields  indicates  that  there  are  velocity  limits  within  which 
scaling  could  be  accomplished.  The  velocity  distributions  around  an  airfoil  will  normally  be 
preserved  up  to  the  stall  as  long  as  the  free-stream  Reynolds  number,  Re,  is  greater  than  2.0 
X  105  (Ref.  1).  This  corresponds  to  a  velocity  of  31.4  fps  for  a  model  airfoil  section  with  a 
chord  of  1.0  ft  and  a  static  temperature  of  0°F  at  an  altitude  pressure  of  10,000  ft.  An  upper 
velocity  limit  is  also  imposed  by  the  onset  of  local  supersonic  flow  on  the  body  at  the  critical 
Mach  number,  Mc.  The  significant  changes  in  the  flow  that  occur  above  this  upper  velocity 
limit  eliminate  the  possibility  of  extrapolating  results  obtained  at  Mach  numbers  below  Mc 
to  above  Mc. 


From  this  discussion,  the  requirements  for  producing  similar  flow  field  characteristics  in 
the  impingement  region  are  postulated  to  be  (1)  the  objects  must  have  scaled  geometries,  and 
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(2)  the  test  Mach  numbers  must  be  greater  than  that  corresponding  to  Re  =  2.0  x  10s  and 
less  than  Mc 

3.2  DROPLET  TRAJECTORY  SCALING  REQUIREMENTS 

Not  only  must  the  flow  about  the  scaled  geometries  be  similar,  but  the  trajectories  of  the 
droplets  entrained  in  this  flow  must  also  be  similar.  A  droplet  scaling  parameter  is  required 
to  produce  geometrically  coincident  regions  of  impingement  and  identical  mass  distributions 
on  the  surfaces  of  two  scaled  bodies.  The  equations  of  particle  motion  must  be  evaluated  to 
identify  the  parameters  that  can  be  used  to  fulfill  these  requirements. 

The  classical  particle  trajectory  equations  can  be  used  to  calculate  droplet  trajectories 
around  an  object  provided  that  an  accurate  flow  field  has  been  previously  calculated.  In 
doing  so,  it  is  assumed  that  the  concentration  of  droplets  is  small,  less  than  700  particles/cc, 
thereby  allowing  two-phase  flow  effects  and  droplet  collision  and  breakup  to  be  neglected. 
The  differential  equation  describing  droplet  trajectories  normally  used  in  icing  studies  is 

(-fr)  -  5  .<» 


A  summary  of  the  simplifications  made  in  the  formulation  of  this  equation  is  given  in  Ref.  2. 


When  an  appropriate  expression  for  the  drag  force,  D,  is  substituted  and  Eq.  (1)  is 
nondimensionalized,  the  following  equation  results: 


(2) 


The  two  nondimensional  terms  are  the  inertia  parameter,  K,  and  the  droplet  Reynolds 
number,  R,  which  can  be  expressed  by  the  following  equations: 


K  = 


_  2  q  w  rj  Uc 


=  gad  1  ua  -  ud  | 

fa 


(3) 

(4) 


The  derivation  of  Eq.  (2)  from  Eq.  (1)  is  shown  in  Appendix  A. 
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Equation  (2)  is  used  to  calculate  the  droplet  trajectory  from  far  upstream  to 
impingement  on  the  body  and,  if  impingement  occurs,  the  total  and  local  collection 
efficiencies.  The  total  collection  efficiency  is  defined  as  the  ratio  of  the  actual  mass  of 
impinging  water  to  the  maximum  value  that  would  occur  if  the  droplets  followed  straight- 
line  trajectories.  Figure  1  illustrates  that  this  definition  can  be  given  in  equation  form  as 

E  =  y° 

-  H  (5) 

where  y0  is  the  vertical  distance  between  the  droplet  release  points  of  the  upper  and  lower 
surface  tangent  trajectories.  The  local  collection  efficiency,  0,  is  also  defined  in  Fig.  1  and 
can  be  written  in  differential  form  as 


0  = 


dyp 

ds 


It  is  related  to  the  total  collection  by  the  equation 

Era  =  lT  I* 


(6) 


(7) 


Su  =  Upper-Surface  Impingement  Limit 
S1  -  Lower-Surface  Impingement  Limit 
H  =  Forward  Projection  of  the  Airfoil  Height 


Total  Collection  Efficiency 

E, 

Local  Collection  Efficiency 

0 


Vo 

H 

ds 


0  ds 


Figure  1,  Definition  of  total  and  local  collection  efficiency. 


Langmuir  and  Blodgett  (Ref.  3)  presented  the  collection  efficiencies  (total  and  local)  and 
the  droplet  impingement  limits  about  simple  geometries  as  functions  of  K  and  R,  as  shown  in 
Figs.  2,  3,  and  4.  They  also  developed  a  term  called  the  modified  inertia  parameter,  Kp,  that 
can  be  expressed  as 
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0.30.01  4  8  30  40  80 


K 

Figure  2.  Total  collection  efficiency  versus  inertia  parameter  for  a  cylinder  (Ref.  3). 


E 

Figure  3.  Stagnation-point  collection  efficiency  versus  inertia  parameter  for  a 
cylinder  (Ref.  3). 
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Ko  =  (X/XJ  K  (8) 

where  X/Xj  is  defined  as  the  range  parameter  which  is  a  function  of  droplet  Reynolds 
number.  A  further  discussion  of  the  droplet  range  parameter  is  given  in  Appendix  A.  When 
the  calculated  impingement  characteristics  are  plotted  as  functions  of  Ko,  the  family  of 
curves  in  Figs.  2,  3,  and  4  collapse  into  approximately  a  single  curve.  Figure  5,  obtained 
from  Ref.  4,  shows  the  plot  of  Em  versus  K0  for  a  cylinder  along  with  the  curves  for  several 
other  geometries. 


Figure  4.  Droplet  impingement  limit  versus  inertia  parameter  for  a  cylinder  (Ref.  3). 

Some  of  the  functional  dependencies  of  Ko  are  illustrated  in  Fig,  6  which  shows  lines  of 
constant  Ko  on  a  plot  of  droplet  size  versus  velocity.  This  plot  was  constructed  by  calculating 
Ko  for  the  specified  full-scale  condition  and  then  calculating  the  droplet  size  required  for  Ko 
to  be  constant  at  various  model  velocities.  The  curves  indicate  that  for  increasing  velocities, 
the  subscale  droplet  size  must  be  reduced  so  that  the  effects  of  droplet  inertia  on  the 
trajectories  are  similar.  Also,  since  a  subscale  model  will  disturb  the  flow  less  than  the  full- 
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scale,  the  droplet  diameter  for  the  model  condition  must  also  be  reduced  to  maintain  the 
effects  of  droplet  inertia.  An  increase  in  static  pressure  is  shown  in  Fig.  6  to  slightly  increase 
the  required  subscale  droplet  diameter.  In  the  range  of  static  temperatures  encountered  in 
most  icing  analysis,  the  effect  of  temperature  on  Ko  and  the  resulting  subscale  droplet 
diameter  is  negligible. 


0.01  0.1  1  10 

Modified  Inertia  Parameter, 


Figure  S.  Total  collection  efficiency  versus  modified  inertia  parameter 
for  various  geometries  (Ref.  4). 

To  verify  Ko  as  a  droplet  scaling  parameter,  it  must  be  shown  that  (1)  the  trajectory 
equations  from  which  Ko  is  derived  can  be  used  to  predict  experimentally  obtained 
impingement  characteristics,  and  (2)  these  characteristics  are  indeed  scaled  for  conditions 
with  constant  Ko-  The  first  requirement  has  been  evaluated  by  previous  droplet  trajectory 
studies  conducted  by  Bragg  (Ref.  2)  and  Frost  (Ref.  5).  Good  agreement  with  experimental 
data  was  obtained  in  both  studies.  Bragg  also  satisfied  the  second  requirement  by  showing 
examples  of  similar  droplet  impingement  characteristics  obtained  at  conditions  with 
constant  values  of  Kg. 

In  a  test  facility,  as  in  nature,  an  icing  cloud  is  composed  of  a  distribution  of  droplet 
diameters.  T o  avoid  calculating  the  value  of  Kq  and  local  collection  efficiency  for  each  of  the 
droplet  sizes,  the  icing  cloud  is  characterized  by  a  single  mass  median  droplet  diameter.  The 
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Full-Scale  Condition:  Kq  =  0.55 
2.0-in.-diam  Cylinder 
Ps  =  12.2  psia 
Ts  =  20°F 
dm  =  20  ym 


Velocity,  fps 

Figure  6.  Effect  of  velocity,  static  pressure,  and  model  size  on  the 
mass  median  droplet  diameter  for  a  constant  value  of 
modified  inertia  parameter,  Kg. 
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mass  median  droplet  diameter  is  defined  as  that  diameter  below  which  half  of  the  total  mass 
of  liquid  water  is  contained.  The  validity  of  this  characterization  has  been  demonstrated  by 
Chang  (Ref.  6).  Figure  7  shows  a  comparison  of  the  local  collection  efficiency  calculated  for 
a  multidispersed  droplet  distribution  (Langmuir  D)  and  for  the  mass  median  droplet 
diameter  characterizing  that  distribution.  Note  that  excellent  comparison  with  experiment  is 
achieved  using  either  the  multidispersed  or  monodispersed  clouds  in  the  stagnation  region 
but  that  the  entire  droplet  distribution  must  be  included  to  accurately  evaluate  the 
characteristics  near  the  impingement  limits.  As  will  be  discussed  in  Section  3.4,  the 
agreement  in  the  stagnation  region  is  sufficient  to  verify  the  use  of  the  mass  median  droplet 
diameter  to  characterize  the  icing  cloud  when  applied  to  scaling  studies.  All  further 
references  to  droplet  diameter  contained  in  the  text  imply  the  mass  median  droplet  diameter 
of  the  icing  cloud  unless  otherwise  specified. 

Sym 

▲  Experimental  Data 
(NACA  TN  3839) 

—  — —  Monodispersed  Droplet 
(dm  =  16.7  pm ) 


Upper  Surface  Lower  Surface 

Distance  Along  Airfoil  Surface  s/c 


Figure  7.  Comparison  of  local  collection  efficiency  calculations  using  monodispersed  and 
multidispersed  droplet  distributions  (Ref.  6). 
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The  above  results  show  that  Kg,  based  on  the  mass  median  droplet  diameter,  satisfies  the 
requirements  of  a  droplet  scaling  parameter.  The  scaling  equation  for  droplet  trajectories 
can,  therefore,  be  written  as  follows: 


(Ko)mode]  —  (Ko)  full-scale 


3J  TOTAL  WATER-CATCH  SCALING  REQUIREMENTS 


(9) 


Droplet  trajectory  scaling  has  assured  that  the  normalized  distribution  of  impinging  mass 
on  the  surface  of  the  body  is  similar.  For  ice  accretions  formed  on  two  geometrically  scaled 
objects  to  be  similar,  the  total  mass  of  impinging  water  per  unit  area  should  be  scaled,  i.e. 


(mw  )model —  h  (mw)  full -scale 


where  k  is  the  model  scale  factor  given  by  the  following  ratio: 

cm 


k  = 


Cf 


(10) 


(11) 


The  full-  and  subscale  characteristic  lengths,  cm  and  cr,  respectively,  can  be  any  scaled 
dimension  but  were  taken  to  be  the  airfoil  chord  or  cylinder  diameter  in  this  study. 
Substituting  Eq.  (11)  into  Eq.  (10),  the  following  equation  results: 


/  mw  \  _  /  \ 

\  c  /model  \  C  / 


full-scale 


(12) 


The  total  mass  of  water  impinging  at  a  specific  location  on  the  surface  in  a  time  interval, 
r,  is  given  by  the  following  equation: 


m;  =  LWC  (U„)  0  t 


(13) 


Substituting  Eq.  (13)  into  Eq.  (12)  results  in  a  relationship  that  will  scale  the  total  mass  of 
water  per  unit  area  impinging  at  each  point  on  geometrically  scaled  objects, 


LWC  (U»)  0  t  _  LWC  (Up,)  0  r 

t  model  C 


full-scale 


(14) 


Since  the  modified  inertia  parameter  scaling  has  ensured  that  the  local  collection  efficiencies 
at  each  location  on  the  surface  are  the  same,  0  can  be  cancelled  from  both  sides  of  Eq.  (14) 
to  yield  the  following  expression: 
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LWC  (U.)  r  =  LWC  (U»)  r 

c  I  model  I  C 

Note  that  Eq.  (15)  is  now  independent  of  the  location,  s,  on  the  surface. 

Past  icing  analyses  by  various  authors  have  related  the  terms  of  Eq.  (15)  to  an  ice  thickness 
by  dividing  by  the  ice  density.  The  resulting  term,  called  the  accumulation  parameter,  can 
therefore  be  expressed  as 


full-scale 


(15) 


Ac  =  .LWCjUoo)  r 
6i  c 

Since  the  type  of  ice  formed  on  a  subscale  accretion  is  to  be  the  same  as  that  formed  on  the 
full-scale,  the  ice  densities  must  be  identical.  Equation  (15)  can  therefore  be  rewritten  in 
terms  of  the  accumulation  parameter  as 


(Ac)model  —  (Ac  )fuii-$cale  (17) 

Figure  8  shows  the  functional  dependencies  of  the  terms  in  the  accumulation  parameter. 
The  plot  shows  that  for  a  constant  LWC,  as  the  velocity  is  increased,  the  icing  time  must 
decrease  so  that  the  total  mass  of  impinging  water  remains  scaled.  Reducing  the  model  size 
at  a  constant  velocity  is  also  seen  to  reduce  the  icing  time  so  that  a  scaled  mass  of  water  is 
deposited. 

In  summary,  this  analysis  has  indicated  that,  if  the  accumulation  parameter,  A<.,  is  held 
constant,  the  total  mass  of  water  per  unit  area  impinging  on  the  surface  of  full-  and  subscale 
geometries  will  be  scaled.  The  scaling  equation  expressing  this  requirement  is  given  by  Eq. 
(17). 

3.4  THERMODYNAMIC  SCALING  REQUIREMENTS 

Thermodynamic  scaling,  or  energy  equation  scaling  as  it  is  commonly  called,  is  required 
so  that  the  type,  surface  characteristics,  and  density  of  the  ice  accretion  is  similar  between 
the  full-  and  subscale  test  articles.  To  formulate  the  mathematical  expression  required  to 
define  the  thermodynamic  scaling  parameters,  it  is  necessary  to  develop  a  physical  model  of 
the  ice-accretion  process. 

The  assumption  made  in  the  discussion  of  the  flow  field  requirements  was  that  if  the 
geometries  are  scaled,  the  normalized  pressure  and  temperature  distributions  in  the  region  of 
droplet  impingement  will  be  similar.  Also,  the  droplet  impingement  parameter,  Ko,  has  been 
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shown  in  Refs.  2  and  5  to  accurately  duplicate  the  distribution  of  mass  on  the  surface  for  two 
scaled  conditions.  Since  these  parameters  serve  as  the  input  to  the  thermodynamic  analysis, 
the  distribution  of  the  thermodynamic  characteristics,  for  example,  the  normalized  surface 
temperature  distribution,  should  also  be  similar  to  that  of  the  full-scale.  To  produce 
identical  thermodynamic  characteristics  in  the  impingement  region  given  similar  input 
distributions,  it  is  assumed  to  be  only  necessary  to  match  the  characteristics  at  one  point. 
The  stagnation  point  is  chosen  because  data  exist  for  the  heat-transfer  coefficient  and  local 
collection  efficiency  in  that  region  for  a  variety  of  body  geometries.  This  assumption  could 
not  be  evaluated  a  priori  because  sufficiently  accurate  experimental,  and  analytical  data  does 
not  exist  at  this  time.  Therefore,  as  in  the  analysis  of  the  flow  field,  the  validity  of  the 
assumption  was  determined  by  comparing  the  full-  and  subscale  ice  accretions.  All  values 
of  the  thermodynamic  parameters  reported  herein  are  evaluated  at  the  stagnation  point.  Recall 
that  Fig.  7  showed  that  the  mass  median  droplet  diameter  could  be  used  to  characterize  the 
droplet  distribution  in  an  icing  cloud  in  the  vicinity  of  the  stagnation  point. 

Full-Scale  Condition,  Ac  =  0.9 
2 . 0-in . -diarn  Cylinder 
LWC  =  0.5  g/m^ 


Velocity ,  fps 

Figure  8.  Effect  of  velocity  and  model  size  on  the  icing  time  for 
a  constant  value  of  accumulation  parameter,  Ac. 
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The  physics  of  stagnation-line  icing  has  been  observed  at  NASA  Lewis  Research  Center 
using  microscopic  observation  techniques  (Ref.  7).  The  observations  have  shown  that  when 
droplets  initially  strike  an  ice-free  surface,  they  coalesce  into  larger  drops  and  then  run  back 
along  the  surface.  By  applying  the  premise  in  Section  2.1  which  states  that  if  the  ice 
accretions  can  be  started  under  similar  conditions,  they  will  grow  similarly,  the  observed 
initial  phase  of  the  ice-accretion  process  is  all  that  has  to  be  modeled,  and  the  more 
complicated  problem  of  modeling  the  formation  of  ice  on  a  previously  deposited  ice  layer  is 
avoided.  A  thermodynamic  equation  of  the  initial  phase  of  the  stagnation-line  process  can 
be  formulated  from  the  above  observations  and  used  to  define  the  scaling  parameters  to  be 
used  in  this  study. 

3.4.1  Icing  Thermodynamic  Model 


The  thermodynamic  analysis  of  an  icing  surface  was  first  presented  by  Tribus  (Ref.  8) 
and  was  used  to  calculate  the  heating  requirements  for  icing  protection  and  proposed  LWC 
measurement  systems.  Messinger  (Ref.  9)  developed  the  model  further  to  include  an  analysis 
of  the  temperature  of  an  unheated  surface  in  icing  conditions  in  three  free-stream 
temperature  regimes,  i.e.  less  than  32°F,  equal  to  32°F,  and  above  32°F,  and  the  concept  of 
the  freezing  fraction,  n,  to  be  discussed  later.  These  early  formulations  have  been  used  in 
various  icing  applications  over  the  years.  Modifications  to  these  early  models  have  consisted 
mainly  of  adding  the  effects  of  compressibility  in  the  flow  about  an  airfoil.  The  model  for 
stagnation-line  icing  used  in  this  study  is  slightly  different  from  that  of  previous  work 
because  both  the  mass  and  energy  balance  are  addressed.  The  form  of  the  energy  equation 
used  in  this  study  can  be  expressed  as  follows: 


The  complete  derivation  of  Eq.  (18)  is  included  in  Appendix  B. 


Dividing  Eq.  (18)  to  the  convective  heat-transfer  coefficient,  hc,  results  in  the  following 
form  of  the  energy  equation: 

-b*  -  B  +  (1  -  n)  -^=-  [Cpw  9ur  (Tsur  -  32)] 
nc 

+  n  [cK  SUT  (Tsur  -  32)  -  Lf]  (19) 
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where 


b  = 


cpws 

he 


(20) 


tf>  =  32  -  Ts  - 


2&  JCpw.s 


0  = 


-T, 


ui,  1  _  we  T 
h,  ^ 


(21) 


(22) 


The  above  terms  will  be  discussed  in  Section  3.4.3.  Before  applying  the  above  equation  to  an 
icing  analysis  it  was  necessary  to  evaluate  experimentally  the  range  of  applicability  of  the 
stagnation-line  thermal  model. 


3.4.2  Verification  of  the  Thermal  Model 


The  experimental  verification  of  the  thermal  model  for  stagnation-line  icing  is  difficult 
and  has  been  limited  to  comparison  between  measured  and  calculated  surface  temperatures, 
icing  onset  temperatures,  and  stagnation-line  ice  growth  rates.  Since  the  experimental 
determination  of  the  icing  onset  temperature  is  somewhat  subjective,  and  the  stagnation-line 
growth  rate  is  dependent  on  precisely  knowing  the  ice  density,  the  comparison  of  surface 
temperatures  is  considered  the  most  accurate  means  of  verifying  the  thermal  model.  Recent 
work  in  this  area  has  been  done  by  Fergus  (Ref.  10),  who  found  that  comparisons  with 
experiment  were  fairly  accurate  when  a  correction  factor  was  applied  to  the  evaporative 
energy  flux  term.  Figure  9  shows  that  the  current  model  compares  well  with  Fergus’ 
corrected  model  without  the  need  for  the  evaporation  term  correction  and  that  both  models 
predict  the  surface  temperature  near  32° F  to  within  1.5°F.  Given  the  errors  in  measuring  the 
temperature  of  a  surface  undergoing  icing,  these  results  provide  a  measure  of  confidence  in 
the  thermal  model  but  definitely  do  not  provide  the  desired  verification.  It  will  therefore  be 
assumed  that  the  accuracy  of  the  thermal  model  is  sufficient  for  application  in  a  scaling 
analysis  and  that  the  validity  of  this  assumption  will  be  determined  by  comparing  the  full- 
and  subscale  ice  accretions. 

3.4.3  Thermodynamic  Scaling  Parameters 

Several  thermodynamic  scaling  parmeters  have  been  proposed  in  past  icing  scaling 
studies  but,  in  general,  without  discussion  of  their  physical  significance.  The  following 
sections  will  define  the  scaling  parameters  used  in  this  study  and  discuss  their  significance 
through  the  presentation  of  ice-accretion  examples  and  figures. 
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Calculated  Surface  Temperature,  °F 

Figure  9.  Comparison  of  measured  and  calculated  surface  temperatures  in 
incipient  idng  conditions. 

3. 4.3.1  Freezing  Fraction 

The  freezing  fraction,  n,  was  defined  by  Messinger  as  being  the  fraction  of  impinging 
liquid  which  freezes  in  the  region  of  impingement.  For  colder  icing  conditions,  the  droplets 
tend  to  freeze  immediately  on  impact  resulting  in  a  rime  ice  accretion,  as  shown  in  Fig.  10a. 
Rime  ice  is  milky  white,  opaque,  and  characterized  by  a  calculated  stagnation-point  freezing 
fraction  equal  to  1.0.  Since  the  droplets  freeze  on  impact,  the  accretion  is  fairly  smooth  and 
aerodynamically  shaped. 

Freezing  fractions  close  to  zero  characterize  glaze  ice  accretions  which  are  composed  of 
clear  ice  and,  in  general,  have  well-defined  horns,  as  shown  in  Fig.  10b.  These  horns  form 
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initially  because  of  Liquid  water  running  along  the  surface  away  from  the  stagnation  point 
and  then  freezing  at  some  point  further  aft  where  the  convective  heat  transfer  is  sufficiently 
increased  to  remove  the  latent  heat  of  fusion.  The  roughness  of  this  accretion  then  enhances 
the  convective  heat  transfer,  causing  the  horns  to  continue  to  grow.  As  the  horns  grow,  the 
collection  efficiency  in  this  region  also  increases,  contributing  further  to  the  growth  of  the 
horns. 


Glaze 

Horns 


b.  Glaze  ice  accretion,  n  -  0.15 


Icing  Conditions 

a 

_b_ 

c 

Cylinder  Diam,  in. 

1.0 

2.0 

1.0 

Rime  Feathers— \ 

Pa.  Psia 

14.2 

12.2 

14.2 

Ts.  °V 

Uc,  fpa 

3m,  pm 

-5.0 

200 

20 

21.1 

367 

22 

5 

400 

15 

Glaze  m 

Center  ll  J 

LWC,  g/m3 

0.6 

0.35 

0.53 

t  ,  min 

8.0 

30.6 

6.  B 

c.  Mixed  ice  accretion,  n  =  0.53 

Figure  10.  Examples  of  ice  accretions  formed  at  various  values  of  stagnation-point 

freezing  fraction. 


Ice  accretions  with  calculated  freezing  fractions  between  approximately  0.3  and  1.0  will 
normally  have  some  of  the  characteristics  of  both  rime  and  glaze  ice  and  are  referred  to  as 
mixed  ice  accretions.  Fig.  10c  shows  that  the  center  portion  will  be  predominantly  e1a7P  but 
will  be  surrounded  by  rime  accretions.  The  size  of  the  glaze  center  will  decrease  as  the 
freezing  fraction  nears  1.0. 

The  predominant  functional  dependencies  of  the  freezing  fraction  are  shown  by  the  lines 
of  constant  freezing  fraction  in  Fig.  11.  As  the  LWC  is  increased,  the  energy  added  to  the 
surface  by  the  droplets  is  also  increased.  The  velocity  must,  therefore,  be  reduced  to 
decrease  the  aerodynamic  heating  and  kinetic  energy  of  the  droplets  to  allow  the  same  fraction 
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of  the  impinging  water  to  freeze.  An  increase  in  the  static  temperature  at  a  constant  velocity 
requires  that  the  LWC  be  reduced,  again  to  maintain  the  proper  freezing  fraction.  A  slight 
variation  of  LWC  with  static  pressure  is  also  shown. 

Full-Scale  Condition:  n  =  0.27 
2.0-in. -diam  Cylinder 
Ps  =  12.2  psia 
Ts  =  20°F 
Uoo  =  250  fps 


Figure  11.  Effect  of  velocity,  static  pressure,  and  static  temperature  on  LWC 
for  a  constant  value  of  freezing  fraction,  n. 

The  above  discussion  indicates  that  the  freezing  fraction  does  have  a  distinct  physical 
significance  in  the  formation  of  ice  accretions.  It  has  been  postulated  in  previous  scaling 
studies  that  the  freezing  fraction  must  be  held  constant  to  produce  similar  ice  accretions  and 
has,  therefore,  been  included  in  this  study.  A  third  scaling  equation  can,  therefore,  be 
written  as  follows: 


(n)  model  —  (^)  full-scale 
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3 .4.3 .2  Relative  Heal  Factor 


The  relative  heat  factor,  b,  was  defined  by  Tribus  (Ref.  8)  as  being  a  measure  of  the  ratio 
of  the  sensible  heat-absorbing  capacity  of  the  impinging  water  per  unit  of  surface  area  to  the 
unit  convective  heat-dissipating  capacity  of  the  same  surface.  It  was  previously  defined  in 
Eq.  (20),  rewritten  here  as 

LWC  (U«)  P  c^ , 


where  c^  s  is  the  specific  heat  of  water.  Note  that  the  relative  heat  factor  is  nondimensional, 
and  since  it  contains  the  local  convective  heat-transfer  coefficient,  1^,  and  local  collection 
efficiency,  0,  is  dependent  upon  geometry. 

Figure  12  shows  lines  of  constant  b  for  the  specified  conditions.  As  indicated  by  Eq.  (24), 
as  the  velocity  increases,  the  LWC  must  decrease.  The  figure  also  shows  that  b  is  relatively 

Full-Scale  Condition,  b  =  0.52 
2. 0-in. -diam  Cylinder 
Pfi  =  12.2  psia 
Ts  =  20°F 
Uoo  =  250  fps 


Velocity,  fps 

Figure  12.  Effect  of  velocity,  static  pressure,  and  model  size  on  LWC  for  a  constant 
value  of  relative  heat  factor,  b. 
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insensitive  to  static  pressure  which  enters  only  through  the  convective  heat  transfer  and  local 
collection  efficiency.  A  decrease  in  the  size  of  the  geometry  will  increase  the  convective  heat- 
transfer  coefficient  requiring  that  the  LWC  also  be  increased  to  maintain  a  constant  value 
of  b.  The  effect  of  static  temperature  on  the  LWC  required  for  the  relative  heat  factor  to 
be  constant  is  negligible,  at  least  in  the  range  of  temperatures  applicable  to  icing  studies. 

It  has  been  postulated  in  previous  icing  studies,  Refs.  11,  12,  and  13  that  the  relative  heat 
factor  must  be  held  constant  to  produce  scaled  ice  accretions.  A  fourth  scaling  equation  to 
be  evaluated  in  this  study  can,  therefore,  be  expressed  as  follows: 


(b)niodel  —  (t)  full-scale 


(25) 


3.4.3J  Droplet  and  Air  Energy  Transfer  Driving  Potentials 

Equation  (19)  identified  two  additional  parameters,  0  and  9,  given  by  Eqs.  (21)  and  (22), 
respectively,  which  are  rewritten  below 


Equations  26  and  27  show  that  0  and  9  have  the  units  of  temperature  and  contain  terms 
relating  to  the  droplet  and  air  energy  transfer,  respectively.  Insight  into  their  signifigance 
can  be  gained  by  examining  the  general  equation  for  the  convective  heat  flux  given  as 

qc  =  hc  AT  (28) 

In  this  equation,  AT  is  the  temperature  difference  that  acts  as  the  driving  potential  for  this 
mode  of  heat  transfer.  Since  0  and  6  also  have  units  of  temperature  and  contain  terms 
relating  to  the  energy  transfer  caused  by  the  droplets  and  the  airflow,  they  can,  by  analogy, 
be  considered  to  be  droplet  and  air  energy  transfer  driving  potentials,  respectively. 

Figure  13  shows  the  relationship  between  the  static  temperature  and  the  velocity  that  is 
indicated  by  Eq.  (26).  This  curve  was  constructed  so  that  the  value  of  0  was  constant.  It  is 
interesting  to  note  that  the  total  or  stagnation  enthalpy,  iT,  is  also  constant  along  this  line. 
Therefore,  0  is  a  measure  of  the  total  enthalpy  of  the  impinging  water. 


27 


AEDC-TR-8  5-30 


Equation  (27)  shows  that  9  is  a  function  of  static  pressure,  temperature,  and  velocity. 
Lines  of  constant  6  for  various  temepratures  are  shown  in  Fig.  14.  Note  that  a  5.0-deg 
increase  in  static  temperature  causes  the  static  pressure  to  be  significantly  reduced.  This 
indicates  that  if  0  is  held  constant,  small  changes  in  the  static  temperature  can  require 
significant  changes  in  static  pressure.  Also,  note  that  static  pressures  above  standard  sea- 
level  pressure  of  14.7  psia  can  be  calculated. 


Figure  13.  Effect  of  velocity  on  static  temperature  for  a  constant 
value  of  droplet  energy  potential,  0, 

In  summary,  the  following  scaling  equations  will  be  evaluated  in  the  experimental 
portion  of  this  study: 


(IQOmodel  —  (Ko)full-scate 
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(Ac)model  =  (^c)  full-scale 
(n)model  ~  (n)full-scalc 
(b)model  =  (b)full-scak 
(^)modcl  =  (<^)  full-scale 
(^)model  =  (®)full-scalc 


Full-Scale  Condition,  0  =  14.40°F 
Ts  =  20°F 
Ps  =  12.2  psia 
Uro  =  250  fps 

■ . .  ■■  ■■  ■■  Effect  of  U«  on  Pg 

-  Effect  of  U®  and  T.  on  P  : 

Ts  -  25°F  8 


Figure  14.  Effect  of  velocity  and  static  temperature  on  static  pressure  for  a  constant 
value  of  air  energy  potential,  0. 
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4.0  SCALING  METHODS  INVESTIGATED 
4.1  PAST  SCALING  INVESTIGATIONS 

In  Section  3.0,  a  proposed  set  of  icing  scaling  parameters  were  presented.  Various 
combinations  of  the  parameters  have  been  postulated  to  produce  similar  ice  accretions  in 
previous  studies.  Table  1  shows  five  of  the  most  well-known  formulations  of  scaling 
equations  and  the  combination  of  parameters  proposed  by  each.  The  equations  proposed  by 
Armand  (Ref.  13)  are  similar  to  those  derived  in  this  report  and,  of  the  studies  listed,  was  the 
only  one  to  provide  experimental  verification.  However,  the  limitations  of  the  test  facility 
used  in  that  study  did  not  allow  adequate  verification  of  the  parameters,  especially  those 
satisfying  thermal  similitude.  Therefore,  experimental  efforts  were  begun  at  the  Arnold 
Engineering  Development  Center  (AEDC)  and  NASA  Lewis  Research  Center  with  the 
purpose  of  comparing  various  scaling  methods  and  experimentally  verifying  a  set  of  scaling 
equations.  The  remainder  of  this  section  will  describe  the  scaling  methods  evaluated  at  the 
AEDC. 


Table  1.  Summary  of  Past  Scaling  Investigations 


Scaling  Analysis 

Scaling  Parameters 

Ko 

Ac 

n 

b 

* 

9 

Douglas  Aircraft 

Co.,  1954  (Ref.  14) 

X 

X 

Lockheed  Aircraft 

Corp.,  1955  (Ref.  11) 

X 

X 

X 

X 

Boeing  Airplane 

Co.,  1962  (Ref.  15) 

X 

X 

British  Aircraft 

Corp-,  1967  (Ref.  12) 

X 

X 

X 

X 

ONER  A  Modane  Centre, 
France,  1977  (Ref.  13) 

X 

X 

X 

X 

X 

X 

4.2  IDENTIFICATION  OF  THE  SCALING  METHODS 

The  purpose  of  the  experimental  portion  of  this  study  was  to  determine  which  of  the 
proposed  scaling  parameters  are  necessary  and  sufficient  to  produce  scaled  ice  shapes.  All 
scaling  studies  agree  that  the  droplets  must  impinge  on  the  surface  in  geometrically  similar 
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locations  and  that  scaled  masses  of  water  per  unit  area  must  impinge  on  the  surface  to 
produce  scaled  ice  shapes.  Therefore,  Ko  and  A^  are  generally  accepted  as  the  droplet  and 
accumulation  scaling  parameters,  respectively.  The  studies  differ  on  which,  if  any,  of  the 
energy  equation  scaling  parameters  are  necessary  to  produce  similar  ice  accretions.  Four 
scaling  methods  were  selected  for  evaluation  in  this  study  and  are  described  below. 

4.2.1  Method  1:  K„  and  At  Constant 

This  method  is  the  simplest  to  apply  because  only  K0  and  A^  are  held  constant.  The  static 
pressure  and  temperature  are  held  constant,  and  the  LWC,  icing  time,  and  velocity  can  be 
varied  as  long  as  Ac  remains  constant.  The  droplet  diameter  is  chosen  so  that  Ko  is  constant. 

4.2.2  Method  2:  Ko,  Ac,  and  n  Constant 

Energy  equation  scaling  is  applied  in  this  method  by  requiring  that  the  freezing  fraction 
be  held  constant  in  addition  to  the  parameters  in  Method  1 .  The  static  pressure,  velocity, 
and  LWC  of  the  model  condition  are  specified,  and  the  static  temperature  is  then  calculated 
so  that  n  is  held  constant. 

4.2.3  Method  3:  K«,  Ac,  n,  and  b  Constant 

In  this  method,  the  relative  heat  factor  is  held  constant  in  addition  to  those  in  Method  2. 
The  velocity  and  static  pressure  of  the  model  condition  can  be  specified,  and  the  LWC  is 
calculated  so  that  b  is  held  constant.  The  static  temperature  of  the  model  icing  condition  is 
again  calculated  from  the  energy  equation  so  that  n  is  held  constant.  The  value  of  Ac  is 
maintained  by  selecting  the  proper  icing  time. 

4.2.4  Method  4:  Ko,  Ae,  n,  0,  and  0  Constant 

The  air  and  droplet  energy  transfer  driving  potential  terms  are  held  constant  in  this 
method  along  with  those  in  Method  2.  It  is  the  most  restrictive  of  the  four  because  only  the 
model  scale  and  velocity  can  be  specified.  All  other  test  parameters  are  determined  from  the 
equations.  A  summary  of  the  methods  and  the  requirements  of  each  are  given  in  Table  2. 

4.3  SOLUTION  OF  THE  SCALING  EQUATIONS 

Many  of  the  past  scaling  investigations  solved  the  scaling  equations  graphically.  The 
complexity  of  the  scaling  equations  contained  in  the  above  methods  make  them  suited  for 
solution  by  computer.  An  icing  similitude  code,  SIMICE,  was  written  for  the  Air  Force  by 
Professor  W.  W.  Bowden  of  Rose-Hulman  Institute  of  Technology  in  1981  under  contract 


31 


AEDC-TR-85-30 


with  the  AEDC.  While  the  form  of  the  code  was  not  applicable  to  this  study,  it  did  contain 
subroutines  that  calculated  many  of  the  physical  constants  that  are  required  for  the  solution 
of  Eq.  (18).  The  code  was  modified  so  that  each  of  the  above  scaling  methods  could  be 
solved  individually  and  was  used  to  calculate  test  conditions  in  the  experimental  portion  of 
this  study. 


Table  2.  Summary  of  the  Scaling  Methods  Investigated  at  the  AEDC 


Method 

Parameters  Held  Constant 

Modified 

Inertia 

Parameter 

Accumulation 

Parameter 

Freezing 

Fraction 

Relative 

Heat 

Factor 

Energy 

Terms 

1 

X 

X 

2 

X 

X 

X 

3 

X 

X 

X 

X 

- 

4 

X 

X 

X 

X 

5.0  EXPERIMENTAL  VERIFICATION 

As  previously  stated,  most  of  the  past  scaling  investigations  did  not  present  experimental 1 
data  verifying  the  proposed  scaling  method.  When  experimental  data  were  presented,  it  was 
not  sufficient  to  identify  the  limitations  of  the  method  so  that  the  data  could  be  used  with 
confidence  in  icing  tests.  The  purpose  of  the  experimental  study  discussed  below  was  to 
investigate  icing  scaling  to  provide  a  more  complete  understanding  of  the  application  of  the 
various  scaling  methods  and  to  determine  which,  if  any,  of  the  proposed  methods  produced 
scaled  ice  accretions.  The  following  section  discusses  the  procedures  used  to  select  test 
conditions  and  to  measure  and  compare  the  resulting  ice  accretions. 

5.1  SCALING  TEST  PROCEDURES 

The  experimental  study  was  begun  by  identifying  and  testing  a  set  of  icing  conditions  on 
a  full-scale  geometry.  The  profiles  of  the  ice  accretions  obtained  from  these  full-scale  test 
conditions  were  measured  using  calipers  and  chilled  carpenter’s  profile  gauges.  When  the  ice 
accretion  was  large  enough  to  allow  removal  from  the  surface,  it  was  removed  and  dipped  in 
beeswax,  producing  a  mold  from  which  a  plaster  casting  was  made.  The  casting  provided 
another  record  of  the  ice-accretion  profile  and  of  the  surface  roughness  characteristics. 
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The  icing  similitude  computer  code  developed  for  this  study,  SIMICE,  was  then  used  to 
calculate  the  test  conditions  for  the  subscale  geometry.  Conditions  were  also  calculated  for 
the  full-scale  geometry  to  evaluate  test  parameter  scaling  discussed  in  Section  2.2.  The 
calculated  test  conditions  were  then  run,  and  the  ice-accretion  characteristics  were  again 
measured  and  recorded. 

Small  differences  between  the  desired  test  condition  and  the  actual  condition  run  are 
present  in  practically  all  types  of  testing.  The  effect  of  these  differences  on  the  values  of  the 
scaling  parameters  can  be  calculated,  but  the  effect  of  differences  in  the  scaling  parameters 
on  the  similarity  of  the  ice  accretions  was  not  known  prior  to  the  tests.  Therefore,  all  test 
conditions  were  reevaluated  after  the  test  using  SIMICE  to  determine  values  of  the  scaling 
parameters  for  the  actual  condition  run.  The  relative  similarity  of  the  full-  and  subscale  ice 
accretions  was  then  evaluated  by  comparing  the  profiles  of  the  accretions  and  the  size  and 
locations  of  the  glaze  horns,  or,  in  the  case  of  mixed  ice  accretions,  the  glaze  and  rime 
portions.  When  necessary,  the  profiles  were  drawn  to  the  same  scale  to  facilitate  more 
accurate  comparisons.  Since  the  drag  increase  caused  by  the  shape  and  surface  roughness  of 
the  ice  accretions  could  not  be  measured  in  the  test  cell,  comparison  of  the  surface 
characteristics  was  limited  to  comparison  of  the  experimental  observations  and  plaster  casts 
of-the  ice  accretion.  ...  ...... 

The  scaling  methods  to  be  evaluated  can  only  be  expected  to  scale  the  dimensions  of  the 
ice  accretions  to  within  the  repeatability  of  two  accretions  formed  at  identical  conditions. 
The  repeatability  of  the  dimensions  is  dependent  on  the  precision  of  the  test  facility  and  the 
measurement  technique.  Figure  15  shows  the  repeatability  of  two  glaze  accretions  formed  at 
identical  atmospheric  and  meteorological  conditions  to  be  on  the  order  of  10  percent. 
Therefore,  a  scaling  method  may  be  considered  successful  if  the  type  of  ice  formed  is  similar 
in  appearance  and  the  dimensions  are  scaled  to  within  approximately  10  percent. 

The  above  procedure  was  applied  for  each  of  the  scaling  methods  to  be  evaluated.  The 
tests  were  run  on  circular  cylinders  of  1-  and  2-in.-diam  and  full-,  1/3-,  and  J/6-scale  airfoils 
at  zero  angle  of  attack.  The  airfoil  section  used  for  these  tests  was  a  supercritical  airfoil  with 
a  reflexed  trailing  edge,  as  shown  in  Fig.  16.  All  cylinder  and  airfoil  test  articles  were  made 
of  solid  aluminum.  The  following  sections  describe  the  test  facility  and  the  procedures  used 
to  set  the  conditions  during  a  test. 

5,2  FACILITY  DESCRIPTION 

The  experimental  portion  of  this  study  was  conducted  in  the  ETF  Research  Test  Cell 
(R-1D)  at  the  AEDC.  The  icing  research  test  cell  (Fig.  17)  consists  of  a  flow-metering 
venturi,  plenum  chamber,  water  spray  system,  bellmouth,  removable  connecting  ducts,  and 
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a  test  chamber.  Water  droplets  are  sprayed  into  the  primary  air  stream  through  a  single,  two- 
phase  atomizing  spray  nozzle  located  in  the  plenum  chamber,  upstream  of  the  bellmouth. 
The  bellmouth  terminates  in  a  12-in. -diam  duct  that  directs  the  conditioned  air  to  the  test 
article  in  a  3-ft.-diam  test  section.  A  secondary  air  system  supplies  air  to  the  test  section  that 
encapsulates  the  primary  flow  to  prevent  the  recirculation  of  water  droplets  around  the  test 
article.  The  air  is  then  exhausted  to  the  atmosphere  either  directly  or  through  the  ETF 
exhaust  plant. 


2.0-in. -dian  Cylinder 


Condition  A 


Condition  A 


-  Condition  B 


Dimension 

Condition 

Percent 

Error 

A 

B 

Impingement  Limit 

1.94 

1.80 

7.3 

Maximum  Thickness 

1.09 

0.94 

13.8 

Maximum  Width 

3.95 

3.78 

4.3 

Figure  15.  Example  of  test  cell  and  ice-accretion  measurement  repeatability. 


Full-Sc!ale  Chord  =  10.907  in. 

Figure  16.  Supercritical  airfoil  section  used  in  the  icing  scaling  studies. 
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Figure  17.  Icing  research  test  cell. 


This  study  requires  that  both  the  droplet  size  and  LWC  in  the  test  section  be  accurately 
known  and  controlled.  The  methods  used  for  specifying  these  parameters  are  described 
below. 

5.2.1  Droplet  Size  Calibration 

A  two-phase  atomizing  spray  nozzle  is  used  to  produce  a  cloud  of  water  droplets  in  the 
test  cell.  A  cloud  with  the  desired  mass  median  droplet  size  and  spray  water  flow  rate  is 
obtained  by  setting  the  air  and  water  pressures  supplied  to  the  spray  nozzle.  Each  nozzle 
must,  therefore,  be  calibrated  prior  to  use  in  an  icing  test  so  that  the  droplet  size 
corresponding  to  a  given  set  of  air  and  water  pressures  is  known.  The  droplet  size 
distribution  in  the  cloud  was  measured  using  a  fiber-optics  particle-sizing  system  (FOS) 
developed  at  the  AEDC.  The  FOS  is  an  imaging  device  that  uses  a  laser  beam  as  a  light 
source  and  an  optical  system  to  define  a  probe  volume  in  the  droplet  flow  field  (Fig.  18).  The 
probe  volume  is  focused  onto  a  linear  array  of  sensor  modules.  As  a  droplet  passes  through 
the  probe  volume,  its  shadow  occludes  a  number  of  the  sensors.  The  number  of  sensors 
occluded  is  proportional  to  the  droplet  diameter.  A  thorough  discussion  of  the  operational 
theory  of  the  FOS  is  given  in  Ref.  16. 

The  FOS  has  been  compared  with  other  particle  diagnostic  systems  including  a 
holographic  imaging  system  (Ref.  17).  Figure  19  shows  a  comparison  of  the  mass  median 
diameters  of  droplet  size  distributions  obtained  from  the  FOS  and  the  spray  nozzle 
calibration  curve  obtained  from  holographic  data.  Based  on  this  and  similar  comparisons, 
the  FOS  is  considered  to  provide  accurate  calibrations  of  droplet  distributions  produced  by 
water  spray  nozzles. 
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Figure  19.  Comparison  of  FOS  data  with  spray  nozzle  calibration  (Kef.  17). 
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5.2.2  Determination  of  LWC 

As  with  the  droplet  size  measurement,  the  LWC  is  not  measured  on-line  but  is  set  by 
introducing  the  proper  water  flow  into  the  air  stream  through  the  spray  nozzle.  Since  the  air 
received  from  the  ETF  air  supply  plant  is  essentially  dry,  considerable  evaporation  of  the 
water  could  occur,  thereby  reducing  the  LWC.  A  computer  code  was  developed  at  AEDC  to 
calculate  the  amount  of  evaporation  from  droplets  in  a  spray  (Ref.  18)  and  was  used  to 
determine  the  additional  amount  of  water  that  must  be  added  to  produce  the  required  LWC 
in  the  test  section.  A  comparison  of  methods  used  to  determine  the  LWC  in  a  test  cell  was 
done  by  Stallabrass  (Ref.  19)  and  showed  that  the  AEDC  method  provided  the  best 
comparison  to  measured  values  of  LWC.  Calibration  of  the  icing  research  test  cell  has  been 
achieved  by  comparing  the  LWC  predicted  by  the  AEDC  method  with  the  LWC  calculated 
from  holographic  droplet-size  measurements.  These  results,  illustrated  in  Fig.  20,  show  that 
the  input  LWC  agrees  well  with  the  measured  values. 


Input  LWC,  g/m^ 

Figure  20.  Comparison  of  input  and  measured  LWC  in  the  research  test  cell  (Kef.  17). 

5.3  TEST  PROCEDURES 

As  discussed  in  Section  5.1,  the  aerodynamic  and  meteorological  conditions  to  be  run 
during  a  scaling  test  were  determined  using  the  computer  code,  SIMICE,  so  that  the 
equations  of  the  scaling  method  to  be  evaluated  were  satisfied.  The  spray  nozzle  pressures 
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required  to  produce  the  desired  water  flow  rate  (LWC)  and  mass  median  droplet  diameter 
were  determined  from  the  appropriate  nozzle  calibration  curves. 

After  a  pretest  calibration  of  all  instrumentation  was  performed,  the  test  cell  pressure 
was  reduced  to  the  desired  pressure  altitude  and  primary  and  secondary  inlet  air  was 
admitted  to  the  test  cell  at  the  required  pressure  and  temperature.  Once  the  test  cell  flow 
conditions  had  stabilized  to  steady  state,  the  desired  water  and  air  flows  through  the  spray 
nozzle  were  set.  Approximately  20  sec  were  required  to  stabilize  the  cloud  on  a  desired 
condition. 

After  exposing  the  test  article  to  the  icing  cloud  for  the  desired  length  of  time,  the  spray 
nozzle  and  primary  air  flow  were  stopped  and  the  test  cell  brought  to  atmospheric  pressure. 
The  connecting  duct  section  immediately  in  front  of  the  test  section  (Fig.  17)  was  removed 
and  the  measurements  of  the  ice  accretions  described  in  Section  5.1  were  obtained. 
Approximately  5  min  were  required  for  this  shutdown  and  measurement  procedure.  Since 
the  test  cell  remained  cold,  no  significant  melting  of  the  ice  accretion  was  observed.  This 
procedure  was  repeated  for  each  test  condition. 

5.4  TEST  RESULTS 

The  scaling  methods  previously  discussed  were  evaluated  beginning  with  the  least 
restrictive,  Method  1 ,  and  progressing  to  the  most  restrictive.  Method  4.  The  ice-accretion 
profiles  included  here  were  redrawn  to  the  same  scale  to  allow  direct  comparisons.  The 
dimensions  shown  on  the  profiles  are  those  that  were  measured  on  the  ice  accretion. 

5.4.1  Method  1:  Kq  and  Ac  Constant 

The  static  pressure  and  temperature  were  held  constant  for  these  cases,  and  the  LWC  and 
icing  time  were  varied  so  that  the  value  of  A<.  was  maintained.  The  droplet  diameter  was 
calculated  so  that  K0  was  constant.  Figure  21a  compares  rime  ice  accretions  formed  on  a 
2-in.-diam  full-scale  cylinder  and  a  l-in.-diam  model  at  a  static  temperature  =  0.0°F.  The 
freezing  fractions  for  these  two  cases  are  approximately  1.0  which  is  indicative  of  rime  ice 
accretions.  The  ice  shapes  are  similar  indicating  that  this  method  appears  to  adequately  scale 
results  for  rime  ice  accretions. 

The  accretions  shown  in  Fig.  21b  were  also  obtained  at  0.0°F  but  at  a  greater  velocity 
than  the  conditions  in  Fig.  21a.  These  accretions  differ  in  shape,  especially  in  the  size  of  the 
rime  feathers  and  glaze  center  portion.  Note  that  the  freezing  fraction  for  the  subscale  model 
is  greater  than  that  of  the  full-scale  model.  This  indicates  that  the  subscale  ice  accretion 
should  have  a  more  rime-like  character  than  the  full-scale.  This  premise  is  evidenced  by  the 
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c.  Velocity  =  200  fps,  mixed  Ice  accretions,  test  article  size  scaling 
Figure  21.  Comparison  of  rime,  mixed,  and  glaze  ice-accretion  scaling  using 
Method  1  to  establish  the  experimental  similitude  conditions. 


larger  rime  accretions  on  both  sides  of  the  glaze  center.  Similar  conclusions  can  be  drawn 
from  the  results  shown  in  Fig.  21c.  Note  that  the  values  of  <j>  and  9  are  the  same  for  the 
comparisons  in  Figs.  21b  and  21c.  This  indicates  that  these  two  parameters  are  not  sufficient 
for  producing  scaled  ice  shapes  but  still  may  be  necessary,  as  will  be  examined  in  the 
evaluation  of  Method  4. 


Figure  21d  shows  glaze  accretions  formed  on  two  full-scale  airfoils.  In  this  set  of 
conditions,  the  velocity  was  scaled  from  400  to  200  fps  and  is  an  example  of  test  parameter 
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scaling.  The  small  amount  of  ice  collected  during  the  400-fps  condition  indicates  that  it  was 
near  the  icing  onset,  as  evidenced  by  the  low  value  of  the  freezing  fraction. 

From  these  results,  it  can  be  concluded  that  adequate  scaling  can  be  achieved  using 
Method  1  only  if  the  conditions  are  cold  enough  to  produce  a  solid  rime  accretion.  If  any 
portion  of  the  accretion  has  a  glaze  character,  this  method  will  not  produce  scaled  ice 
shapes.  These  conditions  also  reemphasized  the  physical  significance  of  the  freezing  fraction 
when  considering  the  type  of  ice  that  will  be  formed. 
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d.  Airfoil  glaze  Ice  accretions,  test  parameter  scaling 
Figure  21.  Concluded. 

5.4.2.  Method  2:  Ko  Ac,  and  n  Constant 

Method  2  differs  from  Method  1  in  that  the  free-stream  static  temperature  of  the  model 
test  condition  is  now  calculated  so  that  the  freezing  fraction  is  held  constant.  Figure  22a 
shows  glaze  ice  accretions  formed  on  two  l-in.-diam  cylinders.  The  LWC  was  scaled  from 
1,2  to  0.8  g/m3  while  holding  the  freezing  fraction  constant.  Both  accretions  are  composed 
entirely  of  glaze  ice,  but  the  size  and  shape  of  the  horns  have  not  been  accurately 
reproduced. 

Figure  22b  shows  mixed  ice  accretions  formed  at  the  same  condition  as  those  in  Fig.  22a 
except  that  the  temperature  was  reduced.  The  profiles  show  some  similarity  at  the  short  and 
intermediate  times,  but  at  the  longest  time,  the  differences  between  the  accretions  have 
become  evident.  The  test  conditions  in  Fig.  22c  also  show  that  the  accretions  appear  similar 
for  short  times  but  become  different  as  the  accretion  grows. 
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These  results  have  shown  that  the  Method  2  can  ensure  that  the  type  of  ice  formed  is  the 
same,  but  the  size  and  locations  of  the  glaze  horns  are  not  adequately  reproduced. 
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a.  Cylinder  glaze  ice  accretions  lest  parameter  scaling 
Figure  22.  Comparison  of  glaze  and  mixed  ice-accretion  scaling  using  Method  2  to  establish 
the  experimental  similitude  conditions. 


5.4.3  Method  3:  K«,  AC1  n,  and  b  Constant 


In  this  method,  the  static  temperature,  droplet  diameter,  and  icing  time  are  determined 
as  in  Methods  1  and  2.  The  LWC  at  which  the  model  is  to  be  tested  is  calculated  so  that  the 
relative  heat  factor  is  constant. 


Figure  23a  shows  two  glaze  accretions  formed  on  a  full-scale  airfoil.  The  static  pressure 
was  held  constant  as  it  would  be  in  an  atmospheric  test  facility,  and  the  velocity  was  scaled 
from  200  to  350  fps.  As  in  Method  2,  the  size  of  the  glaze  horns  is  not  accurately 
reproduced.  Figure  23b  shows  a  condition  for  which  the  static  pressure  of  the  model 
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condition  was  specified  to  be  9.4  psia.  Note  that  this  still  did  not  produce  scaled  ice 
accretions.  Figure  23c,  on  the  other  hand,  shows  mixed  ice  accretions  on  two  l.O-in.-diam 
cylinders  that  were  also  formed  at  conditions  with  different  pressures.  The  scaling  of  these 
two  conditions,  while  not  exact,  is  vastly  improved  over  that  in  Figs.  23a  and  b. 
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b.  Cylinder  mixed  ice  accretions,  test  parameter  scaling 
Figure  22.  Continued. 


These  results  indicate  that,  with  the  proper  selection  of  static  pressure,  this  method  may 
be  able  to  produce  scaled  ice  accretions.  Since  the  additional  knowledge  about  the 
relationship  between  the  velocity  and  static  pressure  is  not  supplied  by  this  method,  accurate 
results  cannot,  in  general,  be  obtained.  It  is  interesting  to  note  that  the  values  of  0  and  $  are 
almost  constant  in  Fig.  23c  indicating  that  these  may  supply  the  additional  information 
required  to  determine  the  correct  static  pressure. 


42 


AEDC-TR-8  5-30 


0.96 


TT^ 

1.08  (  0.82  ] 

iL-xly 


0.91  |*_ 
6.81 


Ac 

0.38 

0. 

76 

1.14 

LWCj 

B/m3 

fps 

Ts. 

°F 

PS. 

psia 

dm, 

ym 

*0 

n 

b 

°F 

1 

CD  O 

A  2.0-in. 

-diam  0.53 

300.7 

-2.3 

14 

26 

2.24 

0.54 

0.98 

32.6 

39.7 

B  1.0-in. 

-diam  0.53 

400.4 

5 

14.2 

15 

2.20 

0.52 

0.51 

23.9 

24.5 

c.  Cylinder  mixed  ice  accretions,  test  article  size  scaling 
Figure  22.  Concluded. 


28  8.5  0.25  0.32  0.13  0.22  fi.O  2.7 


a.  Airfoil  glaze  ice  accretions,  test  parameter  scaling 
Figure  23.  Comparison  of  glaze  and  mixed  ice-accretion  scaling  using  Method  3  to  establish 
the  experimental  similitude  conditions. 
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c.  Cylinder  mixed  Ice  accretions,  test  parameter  scaling 
Figure  23.  Concluded. 
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5.4.4  Method  4:  Ko,  Ac,  n,  <j>,  and  8  Constant 

In  this  method,  the  static  temperature  and  pressure  are  calculated  so  that  $  and  9  are  held 
constant  between  the  model  and  full-scale.  Since  the  other  icing  parameters  are  calculated  as 
in  Method  2,  only  the  model  scale  factor  and  free-stream  velocity  can  be  specified. 


The  results  in  Fig.  24a  show  that  excellent  similitude  has  been  achieved  between  two  glaze 
ice  accretions.  Note  that  the  static  pressures  are  different  between  the  model  and  full-scale 
conditions  and  have  been  calculated  so  that  6  is  held  constant. 
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a.  Cylinder  glaze  ice  accretions,  test  article  size  scaling 
Figure  24.  Comparison  of  glaze  and  mixed  ice-accretion  scaling  using  Method  4  to  establish 

the  experimental  similitude  conditions. 


The  comparison  in  Fig.  24b  illustrates  that  scaled  accretions  can  also  be  obtained  for  test 
parameter  scaling.  Condition  A  was  obtained  at  a  pressure  altitude  of  10,000  ft.  The 
application  of  Method  4  allowed  the  condition  to  be  scaled  down  to  sea  level,  as  shown  in 
Condition  B,  so  that  the  condition  could  be  run  in  an  atmospheric  test  facility. 
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This  method  also  produced  accurate  scaling  results  for  the  mixed  ice  accretions  on  the 
full-  and  1/3-scale  airfoils  shown  in  Fig.  24c.  Note  that  even  though  the  freezing  fraction  is 
calculated  to  be  equal  to  1.0,  the  ice  in  the  center  region  has  some  glaze  characteristics, 
namely  the  well-defined  horns.  This  indicates  that  a  calculated  freezing  fraction  equal  to  1.0 
is  not  sufficient  to  ensure  a  smooth,  aerodynamically  shaped,  rime  ice  accretion. 

Figure  24d  also  shows  well-scaled  results  for  glaze  ice  accretions  formed  on  the  full-  and 
1 /3-scale  airfoils.  The  accretion  on  the  1 /3-scale  model  is  slightly  smaller  because,  upon 
completion  of  the  posttest  analysis,  the  accumulation  parameter  was  found  to  be  less  than 
that  required  for  exact  scaling.  It  is  interesting  to  note  that  the  accumulation  parameters 
differed  by  21  percent,  and  the  maximum  thickness  of  the  ice  accretions  differed  by  15 
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b.  Cylinder  glaze  ice  accretions,  test  parameter  scaling 
Figure  24.  Continued. 
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percent.  This  illustrates  the  importance  of  the  posttest  calculation  of  the  scaling  parameters 
in  the  interpretation  of  scaling  test  data. 


c.  Airfoil  mixed  ice  accretions,  test  article  size  scaling 


d.  Airfoil  glaze  ice  accretions,  test  article  size  scaling 
Figure  24.  Concluded. 


In  all  comparisons  where  this  method  was  applied,  the  type  and  texture  of  the  ice  formed 
was  similar.  Comparison  of  the  dimensions  of  the  scaled  shapes  show  that  most  dimensions 
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are  geometrically  scaled  to  within  10  percent.  As  was  discussed  in  Section  5.1,  since  this 
scaling  method  produces  results  to  the  same  order  of  accuracy  as  the  repeatability  of  the 
measurement  technique,  it  can  be  concluded  that  the  method  accurately  scales  the  flow  field, 
droplet  impingement,  and  thermodynamic  characteristics  of  the  full-scale  icing  condition. 

These  results  have  identified  a  scaling  procedure  that  accurately  duplicates  the  ice- 
accretion  profiles.  In  this  scaling  method,  Ko,  A*.,  n,  0,  and  8  are  held  constant.  The  set  of 
verified  scaling  equations  can,  therefore,  be  expressed  as  follows: 


1-  (f^o)model  —  (Ko)full-scale 

2.  (Ac)mode|  —  (Ac)fu||_scaie 


3*  (n)modcl  —  (n)full-scalc 


4-  (0) model  (0) full-scale 


(^) model  —  (^)full-scale 


A  final  version  of  the  icing  scaling  similitude  code,  S1MICE,  was  written  to  solve  the 
above  set  of  verified  scaling  equations  (See  AEDC-TR-85-30,  Vol.  II).  The  icing  test 
parameters  that  must  be  determined  from  thse  equations  are  as  follows: 

1.  velocity, 

2.  static  pressure, 

3.  static  temperature, 

4.  liquid-water  content, 

5.  droplet  diameter,  and 

6.  icing  time. 

Note  that  there  are  six  icing  test  parameters  and  only  five  equations,  indicating  that  the 
problem  is  underspecified  and  that  one  of  the  test  parameters  can,  therefore,  be  arbitrarily 
selected.  In  practice,  this  additional  parameter  is  selected  to  overcome  a  test  facility 
limitation.  The  solution  of  the  above  equations  will  be  illustrated  by  two  practical  scaling 
examples  to  be  discussed  in  the  following  section. 
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6.0  ICING  SCALING  APPLICATIONS 
6.1  SAMPLE  SCALING  PROBLEMS 

The  experimental  study  has  resulted  in  the  verification  of  Ko,  Ac,  n,  0,  and  6  as  scaling 
parameters.  The  purpose  of  this  section  is  to  graphically  illustrate  the  solution  of  the  scaling 
equations  for  two  icing  test  applications. 

In  these  examples,  it  is  desired  to  simulate  the  ice  accretion  that  would  be  formed  on  a 
36.0-in.  chord  NACA  0012  airfoil  when  exposed  to  the  following  icing  condition: 


geometry 

chord 

leading-edge  diameter 
static  temperature 
static  pressure 
velocity 

droplet  diameter 
liquid-water  content 
icing  time 


=  NACA  0012  airfoil 
=  36.0  in. 

=  1.14  in. 

=  25.0°F 
=  10.1  psia 
=  200.0  fps 

—  20.0  tim 

=  0.5  g/m3 
=  10.0  min 


The  icing  similitude  code,  SIMICE,  was  used  to  evaluate  the  scaling  parameters  for  the 
full-scale  condition  given  above  using  the  equations  presented  in  the  text.  The  results  of 
these  calculations  are  shown  below. 


Ko  =  0.036 
Ac  =  0.013 
n  =  0.260 
0  =  6.190°F 
6  =  9.58 1°F 

Two  hypothetical  test  facilities,  each  of  which  has  a  limitation  that  does  not  allow  the 
full-scale  condition  to  be  run,  will  be  discussed  in  the  following  sample  problems. 

6.1.1  Sample  Problem  1:  Test  Article  Size  Scaling 

The  full-scale  icing  condition  is  to  be  simulated  in  a  test  facility  where,  because  of 
blockage  effects,  the  test  article  can  have  a  maximum  chord  of  12.0  in.  The  scaling  equations 
must  be  applied  to  determine  the  model  icing  conditions  so  that  the  ice-accretion  shape  will 
be  simulated. 
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The  scaling  solution  is  begun  by  producing  the  plot  of  static  temperature  versus  velocity 
shown  in  Fig.  25.  This  curve  is  constructed  from  Eq.  {21)  so  that  $  has  the  same  value  as 
the  full-scale  icing  condition,  namely  6. 19°F.  Since  one  of  the  test  parameters  can  be  specified, 
let  the  velocity  of  the  model  test  condition  be  150.0  fps.  At  this  velocity,  Fig.  25  shows  that 
the  free-stream  static  temperature  of  the  model  condition  should  be  25.4<>F. 

The  static  pressure  is  now  determined  so  that  9  is  constant  between  the  model  and  full- 
scale.  Figure  26  shows  the  static  pressures  and  velocities  for  which  9  is  constant.  Remember 
from  Eq.  (22)  and  Fig.  14  that  9  is  also  a  function  of  the  static  temperature.  Therefore,  the 
static  temperature  is  varying  along  the  curve  in  Fig.  26  as  a  function  of  the  velocity  so  that 
is  constant  and  equal  to  6.19°F  as  shown  in  Fig.  25.  For  a  velocity  of  150.0  fps,  the  static 
pressure  required  for  scaled  ice  shapes  is  found  to  be  12.2  psia.  Note  that  for  velocities  below 
100.0  fps,  static  pressures  greater  than  sea-level  static  pressure  are  calculated. 

The  LWC  can  now  be  determined  so  that  the  freezing  fraction  is  constant.  Figure  27 
shows  values  of  LWC  as  a  function  of  velocity  for  which  the  freezing  fraction  of  the  model 
condition  is  equal  to  that  of  the  full-scale.  Since  both  the  static  temperature  and  pressure  are 


Figure  25.  Sample  Problem  1:  Static  temperature  of  the  model  condition  as  a  function  of 
model  velocity  for  a  constant  value  of  droplet  energy  potential,  <t». 
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Figure  26,  Sample  Problem  1:  Static  pressure  of  the  model  condition  as  a  function  of  model 
velocity  for  a  constant  value  of  air  energy  potential,  0. 

functions  of  velocity,  they  are  also  varying  along  this  curve.  Entering  Fig.  27  with  a  velocity 
of  150.0  fps,  the  LWC  required  for  the  model  test  condition  is  found  to  be  1.10  g/m3. 

Figures  28  and  29  are  used  in  a  similar  manner  to  determine  the  mass  median  droplet 
diameter  and  the  icing  time  so  that  the  modified  inertia  parameter  and  accumulation 
parameter,  respectively,  are  held  constant. 

By  application  of  the  scaling  equations,  the  conditions  at  which  the  model  should  be 
tested  to  obtain  a  scaled  ice  shape  have  been  found  to  be  as  follows. 


velocity 

=  150  fps 

static  pressure 

=  12.2  psia 

static  temperature 

=  25.4°F 

LWC 

=  1.10  g/m3 

droplet  diameter 

=  12.1  (tm 

icing  time 

=  2.02  min 
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If,  upon  completion  of  this  analysis,  it  was  found  that  the  LWC,  say,  was  not  obtainable  at 
the  required  droplet  diameter,  the  procedure  would  be  repeated,  starting  with  Fig.  25,  with  a 
new  choice  for  the  velocity  and  would  be  continued  until  a  set  of  model  conditions  that 
could  be  run  were  obtained. 


Figure  27.  Sample  Problem  1:  LWC  of  the  model  condition  as  a  function  of  model  velocity 
for  a  constant  value  of  freezing  fraction,  n. 

6.1.2  Sample  Problem  2:  Test  Parameter  Scaling 

In  this  example,  the  full-scale  test  condition  is  to  be  simulated  in  an  atmospheric  test 
facility  in  which  the  total  pressure  is  constant  and  equal  to  the  ambient  pressure  of  14.0  psia. 
The  test  facility  is  large  enough  so  that  the  full-scale  geometry  can  be  tested. 
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Since  and  8  are  independent  of  the  model  size.  Figs.  25  and  26  are  also  applicable  to 
this  problem.  Since  there  is  now  a  restriction  on  the  static  pressure,  the  solution  should  be 
started  in  Fig.  26.  In  this  example,  it  would  be  convenient  to  replot  Fig.  26  in  terms  of  the 
total  pressure,  as  shown  in  Fig.  30.  Entering  Fig.  30  at  14  psia  shows  that  a  velocity  of  114 
fps  is  required  for  an  accurate  simulation.  Figure  25  is  then  used  as  before  to  determine  the 
required  static  temperature. 


Velocity,  fps 


Figure  28.  Sample  Problem  1:  Mass  median  droplet  diameter  of  the  model  condition  as  a 
function  of  model  velocity  for  a  constant  value  of  modified  inertia  parameter, 
K0. 

Since  the  velocity,  pressure,  and  temperature  are  different  from  those  in  Sample 
Problem  1,  it  is  necessary  to  calculate  the  LWC,  droplet  diameter,  and  icing  time  as 
functions  of  velocity  to  produce  the  plots  shown  in  Figs.  31,  32,  and  33,  respectively.  These 
are  then  used  in  the  same  manner  as  in  Sample  Problem  1  to  determine  the  remaining  test 
conditions.  The  solution  to  this  sample  scaling  problem  is  as  follows: 
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velocity 
static  pressure 
static  temperature 
LWC 

droplet  diameter 
icing  time 


Figure  29.  Sample  Problem  1:  Icing  time  of  the  model  condition  as  a  function  of  model 
velocity  for  a  constant  value  of  accumulation  parameter,  Ac. 

Since  the  test  facility  has  a  strict  pressure  restriction,  this  is  the  unique  solution  at  the 
required  pressure  that  satisfies  all  the  scaling  equations  for  the  full-scale  article.  If  the 
facility  cannot  obtain  one  of  the  above  conditions,  for  example,  the  droplet  diameter  at  the 
required  LWC,  it  would  be  necessary  to  consider  reducing  the  size  of  the  test  article  and 
attempting  to  find  a  solution  at  the  reduced  size.  It  must  be  realized  that  a  condition  that 
satisfies  the  scaling  equations  and  is  within  the  operational  limits  of  the  test  facility  may  not 
exist. 


=  114  fps 
=  13.9  psia 
=  25.6°F 
=  0.78  g/m* 
=  26.7  jim 
=  1 1.3  min 
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Velocity,  fps 

Figure  31.  Sample  Problem  2:  LWC  of  the  model  condition  as  a  function  of  model  velocity 
for  a  constant  value  of  freezing  fraction,  n. 
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Velocity,  fps 

Figure  32.  Sample  Problem  2:  Mass  median  droplet  diameter  of  the  model  condition  as  a 
function  of  model  velocity  for  a  constant  value  of  modified  inertia  parameter, 
Ko. 


Velocity,  fps 

Figure  33.  Sample  Problem  2:  Icing  time  of  the  model  condition  as  a  function  of  model 
velocity  for  a  constant  value  of  accumulation  parameter,  Ac. 
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The  previous  examples  have  illustrated  the  solution  to  two  practical  scaling  problems. 
The  icing  similitude  computer  code,  SIMICE,  was  used  to  produce  the  figures  shown  in  this 
section. 


7.0  LIMITATIONS 

7.1  VELOCITY  SCALING  LIMITATIONS 

7.1.1  Critical  Mach  Number 


As  discussed  in  Section  3.1,  the  critical  Mach  number  of  the  test  geometry  serves  as  the 
upper  limit  for  velocity  scaling  because  the  character  of  the  flow  field  changes.  This  limit 
also  exists  in  the  droplet  trajectory  calculations  because  the  effects  of  compressibility  are  no 
longer  negligible  (Ref.  20).  Therefore,  scale  accretions  cannot  be  expected  if  the  velocity  is 
scaled  above  Me.  If  velocity  scaling  is  desired  with  both  velocities  greater  than  or  equal  to 
Me,  the  droplet  trajectories  for  each  velocity  should  be  calculated  to  determine  if  there  are 
significant  differences  caused  by  compressibility.  If  necessary,  relationships  that  express  the 
stagnation-line  collection  efficiency  as  a  function  of  modified  inertia  parameter  must  be 
calculated. 

7.1.2  Ice  Shedding 

The  experimental  results  have  shown  that  ice  accretions  with  similar  appearance  can  be 
obtained  through  the  application  of  velocity  scaling.  The  shedding  of  ice  accretions  imposes 
limits  on  velocity  scaling  that  are  more  restrictive  than  the  critical  Mach  number  limit 
described  above. 


Ice  will  shed  from  an  object  when  the  shear  stress  between  the  ice  and  the  surface  or 
within  the  ice  structure  itself  reaches  some  critical  value  which  is  dependent  on  the 
conditions  at  which  the  ice  was  formed.  Since  the  model  and  full-scale  ice  accretions  are 
formed  under  similar  thermodynamic  processes,  the  critical  values  of  the  shear  stress  will  be 
nearly  identical.  The  stress  within  the  ice  is  a  result  of  the  drag  and  moment  forces  on  the  ice 
accretion  itself,  as  shown  in  Fig.  34.  The  magnitude  of  these  forces  are  dependent  on  the  size 
of  the  accretion,  static  pressure,  and  free-stream  velocity. 


Since  the  forces  are  proportional  to  the  dynamic  pressure,  q  =  1/2  y  Ps  M2  =  1/2  ga  U  J, 
an  equivalent  force,  F,  can  be  defined  such  that  the  shear  stress  at  the  surface  of  the  airfoil 
can  be  expressed  as 


a  — 


Ga  U 


2 

00 


Aref 


As 


(29) 
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where  cd;  is  a  total  drag  coefficient  that  can  be  attributed  to  the  ice  accretion  and  is  a 
function  of  the  shape  of  the  ice  accretion.  The  areas  Aref  and  A,  are  defined  in  Fig.  34.  Since 
the  full-  and  subscale  ice  accretions  have  similar  profiles  and  surface  characteristics,  the 
values  of  cDj  will  be  similar  assuming  the  effects  of  the  difference  in  Reynolds  number  to  be 
small. 


C  =  Centroid  of  the  Ice  Accretion  Cross-Sectional  Area 

A  „  =  Reference  Dimension  for  the  drag  coefficient.  Cr, 
ret  D 

A  =  Area  of  Attachment 
s 

F  =  Aerodynamic  Force  Acting  Through  the  Centroid 
M  =  Moment  of  the  Ice  Accretion  About  the  Centroid 

Figure  34.  Definition  of  icing  shedding  analysis  terms. 


Suppose  that  a  full-scale  accretion  has  been  formed  at  a  velocity  U|  and  obtained  a  value 
of  c^  such  that  the  critical  shear  stress  has  just  been  reached,  i.e.  shedding  is  imminent.  The 
shear  stress  at  that  time  is 


6a  LI]  Aref 
As 


(30) 


A  subscale  accretion  will  have  the  same  value  of  Cqj  and  same  ratio  of  Aref/As,  but  since 
it  has  been  formed  at  a  lower  velocity,  the  shear  stress  will  not  be  near  the  critical  value. 
While  many  simplifying  assumptions  have  been  made,  this  analysis  does  show  that  (1) 
testing  at  lower  velocities  can  allow  ice  shapes  to  grow  larger  than  they  would  on  a  full-scale 
geometry,  and  (2)  it  may  be  possible  to  scale  the  shedding  characteristics  by  holding  the 
dynamic  pressure  constant. 


At  the  time  this  experimental  study  was  conducted,  the  icing  research  test  cell  was  not 
capable  of  adequately  evaluating  this  proposed  scaling  requirement  because  of  air  mass  flow 
rate  limitations.  Limited  results  did  reveal  that  above  approximately  M  =  0.4  at  a  static 
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pressure  of  14.2  psia  (q  —  1.6  psia)  shedding  played  a  significant  role  in  the  ice-accretion 
process  and  affected  the  final  shape  of  the  ice  accretion.  Therefore,  without  verification  of 
the  dynamic  pressure  scaling  requirement,  scaling  from  a  full-scale  dynamic  pressure,  qf, 
greater  than  1.6  psia  to  a  model  dynamic  pressure,  qm,  less  than  1.6  psia  or  vice  versa  is  not 
recommended.  Accurate  velocity  scaling  with  both  qf  and  qm  greater  than  or  less  than  1.6 
psia  without  setting  qm  =  qf  can  be  accomplished  only  if  the  time  to  shed  for  either  ice 
accretion  has  not  been  exceeded. 


Analytically  predicting  the  time  to  shed  is  dependent  on  the  ability  to  evaluate  (1)  the 
change  in  drag  and  moment  coefficient  as  functions  of  time,  and  (2)  the  critical  shear  stresses 
between  the  ice  and  the  surface  and  within  the  ice  itself.  These  parameters  are  also  functions 
of  the  aerodynamic  and  meteorological  conditions  at  which  the  accretion  was  formed  and 
are  beyond  the  present  capabilities  of  icing  analysis.  Only  test  results  and  empirical 
correlations  can  currently  provide  direction  in  defining  these  limits. 

7.2  GENERAL  LIMITATIONS 

The  ranges  over  which  the  atmospheric  and  meteorological  parameters  were  tested  are 
shown  in  Table  3.  These  limtis  are  the  result  of  facility  limitations  and  riot  of  the  scaling 


Table  3.  Ranges  of  Test  Parameters  Investigated  at  the  AEDC 


Parameter 

Low 

High 

Velocity,  fps 

100 

400 

Static  Pressure,  psia 

4.4 

14.2 

Static  Temperature,  °F 

-5 

32 

LWC,  g/m3 

0.26 

1.54 

Droplet  Diameter, 

10 

45 

Icing  Time,  min 

1.5 

31.0 

Model  Scale  Factor 

0.17 

6.0 

equations  themselves.  The  test  results  did  not  indicate  that  any  of  these  were  physical  limits 
of  the  scaling  equations.  With  the  possible  exception  of  velocity,  these  ranges  adequately 
cover  a  large  portion  of  the  icing  test  envelopes  required  by  aircraft  manufacturers. 
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8.0  SUMMARY  OF  RESULTS 

The  objectives  of  this  study  were  to  evaluate  the  equations  governing  the  ice-accretion 
process  to  identify  proposed  scaling  parameters  and  to  conduct  tests  to  determine  which,  if 
any,  of  the  proposed  methods  produced  scale  ice  accretions.  A  set  of  scaling  equations  were 
verified,  and  limitations  of  the  seating  techniques  were  identified.  A  computer  code  was 
developed  to  solve  the  scaling  equations  for  use  in  various  icing  applications.  The  results  of 
this  study  can  be  summarized  as  follows: 

1 .  A  set  of  equations  that  can  be  used  to  calculate  test  conditions  so  that  scaled  ice 
shapes  are  produced  on  geometrically  similar  bodies  was  developed  and 
experimentally  verified.  The  test  conditions  are  such  that  K©,  Ac,  n,  <j>,  and  0,  are 
held  constant. 

2.  The  posttest  evaluation  of  the  scaling  parameters  based  on  the  actual  test  conditions 
was  necesary  for  the  accurate  evaluation  of  test  results. 

3.  An  icing  similitude  computer  code,  S1M1CE,  was  developed  to  calculate  similitude 
conditions  using  the  verified  scaling  equations. 

4.  The  equations  are  applicable  over  the  range  of  meteorological  conditions  found  in 
natural  icing,  with  the  possible  exception  of  velocity. 

5.  Velocity  is  the  primary  limitation  of  the  scaling  equations.  To  maintain  scaled  flow 
fields  and  droplet  impingement  characteristics,  both  the  model  and  full-scale 
velocities  must  yield  a  Reynolds  number  greater  than  or  equal  to  2.0  x  105  and  less 
than  the  velocity  giving  the  critical  Mach  number  of  the  body  geometry. 

6.  At  values  of  dynamic  pressure  of  approximately  1.6  psia,  the  shedding 
characteristics  of  the  ice  accretion  were  observed  to  effect  the  final  ice  shape. 
Therefore,  ice  accretions  did  not  scale  well  when  qr  >  1.6  psia  and  qm  <  1.6  or  vice 
versa. 

7.  A  method  for  scaling  shedding  characteristics  by  maintaining  the  dynamic  pressure 
was  developed,  but  facility  limitations  did  not  permit  adequate  verification. 
Additional  testing  is  required  to  verify  this  ice-shedding  scaling  technique. 
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APPENDIX  A 


DERIVATION  OF  THE  MODIFIED  INERTIA  PARAMETER 


The  droplet  trajectory  equation  was  given  in  the  text  as 


(A-l) 


Further  simplication  of  this  equation  is  made  by  assuming  that  the  droplets  are  spherical. 
Approximating  droplets  as  rigid  spheres  is  valid  for  drop  radii  less  than  500  pm  (Ref.  A-l), 

The  drag  force  for  a  sphere  in  Stokes’  flow  (R  <  1.0)  is  given  by  the  equation 


f^Stokes  —  /la  Tj  (Ua  U(j) 


(A-2) 


Since  most  droplet  flow  occurs  outside  this  range,  a  correction  must  be  made  to  account  for 
deviations  from  Stokes*  Law.  This  correction  has  the  form 


D  =  6irMard(Ua-Ud) 


CDR 

24 


(A-3) 


When  Stokes’  Law  is  applicable,  Cp(R)/24  =  1 .0,  and  the  drag  force  is  equal  to  DSu>kes- 
Under  other  conditions,  Cp(R)/24  <  1.0,  and  accounts  for  deviations  from  Stokes’  Law. 


The  effect  of  gravity  on  the  trajectories  of  droplets  less  than  50.0-/tm  is  usually  negligible 
and,  therefore,  is  omitted  in  most  icing  studies  (Ref.  A-2).  With  these  simplications,  the 
droplet  trajectory  equation  becomes 


(A-4) 


Note  that  two  nondimensional  parameters  have  appeared  in  Eq.  (A-4).  The  first  is  the 
droplet  Reynolds  number,  R,  expressed  as 


R  = 


Q&d  I  L'a  —  Ld 
Ma 


(A-5) 


The  second  is  the  inertia  parameter,  K,  which  can  be  expressed  as 

2  Qw  rd  L  oo 


(A-6) 
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where  c  is  the  characteristic  length  of  the  body  geometry,  i.e.  airfoil  chord  or  cylinder  radius. 

Since  all  the  terms  in  the  equation  are  now  nondimensional,  solutions  for  equal  values  of 
K/(CdR/24)  will  be  identical,  assuming  that  the  requirement  for  identical  flow  fields  has 
been  satisfied.  This  term  is  not  constant  throughout  the  trajectory  calculation  since  R  is 
based  on  the  velocity  difference  between  the  gas  and  the  droplet.  Therefore,  a  representative 
value  of  K/(CdR/24)  that  is  unique  for  given  trajectory  will  serve  as  a  droplet  scaling 
parameter. 


Langmuir  (Ref.  A-3)  used  an  average  value  of  CD(R)/24  to  represent  each  trajectory  in 
his  analysis.  This  average  was  defined  as  the  range  parameter,  X/\S1  and  can  be  expressed  by 
the  following  equation: 


X/X,  =  —  f  Ru 

Ru 


24 

C&m 


dm 


(A-7) 


where  Ru  is  the  droplet  Reynolds  number  based  on  the  free-stream  velocity,  U^.  If  values  of 
CD  as  a  function  of  Ru  are  known,  this  equation  can  be  integrated  to  produce  values  of  \/\ 
as  a  function  of  Ru  as  shown  by  the  plot  in  Fig.  A-l .  The  modified  inertia  parameter,  Ko  is 
then  defined  as 


K0  =  (X/Xs)  K  (a-8) 

An  alternate  approach  proposed  by  Bragg  (Ref.  A-2)  used  the  value  of  CDR/24 
corresponding  to  Ru  instead  of  the  average  used  by  Langmuir.  In  this  case,  the  scaling 
parameter  becomes 


K  = 


K 

(CDRu/24) 


(A-9) 


By  substituting  an  expression  for  CD  of  the  form  CD  =  Rj,  Eq.  (A-9)  becomes 

K - ; -  (A-10) 

and  is  solved  by  substituting  the  value  of  7  that  best  fits  the  CD  versus  Ru  curve  in  the  desired 
Reynolds  number  range. 
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Although  both  of  these  methods  have  been  shown  to  adequately  scale  droplet 
trajectories,  the  method  of  Langmuir  was  used  in  the  current  study  because  all  of  the  past 
droplet  impingement  data  have  been  plotted  as  function  Kq. 
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Figure  A-l.  Droplet  range  parameter  versus  Reynolds  number. 
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APPENDIX  B 

DERIVATION  OF  THE  ICING  ENERGY  EQUATION 

DEFINITION  OF  THE  CONTROL  VOLUME 

The  control  volume  to  be  analyzed  is  located  along  the  stagnation  line  and  extends  from 
outside  the  boundary  layer  to  the  surface  of  the  body  as  shown  in  Fig.  B-l.  It  encloses  a 
distance  As  along  the  external  surface  and,  for  dimensional  completeness,  extends  one  unit 
length  in  the  spanwise  direction  (into  the  page). 

MASS  BALANCE  ON  AN  ICING  SURFACE 


An  evaluation  of  all  the  mass  entering  and  Leaving  the  control  volume  is  shown  in  Fig. 
B-2.  A  mass  balance  equation  can  be  formed  from  these  terms  and  is  given  below. 


WwAs  -  WgAs  -  wr  =  WjAs  (B-l) 

There  will  be  no  water  inflow  along  the  surface  because  a  stagnation-point  control  volume 
has  been  selected. 

Since  the  freezing  fraction  is  defined  as  the  proportion  of  the  total  mass  of  liquid 
entering  the  control  volume  that  freezes  in  that  control  volume,  it  can  be  expressed  by  the 
following  equation: 


n  = 


Ww 


(B-2) 


By  substituting  Eq.  (B-l)  into  Eq.  (B-2),  the  water  flow  out  of  the  control  volume  can  be 
expressed  as 


wr  =  (1  -  n)  WwAs  -  WeAs  (B-3) 

ENERGY  BALANCE  ON  AN  ICING  SURFACE 

The  same  control  volume  concept  is  used  to  formulate  the  energy  balance  on  the  icing 
surface.  In  words,  the  First  Law  of  Thermodynamics  for  a  control  volume  can  be  expressed 
as 


energy  inflow  rate  =  energy  outflow  rate  +  energy  storage  rate 
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The  modes  of  energy  transfer,  illustrated  in  Fig.  B-3,  are  as  follows: 


Mode  of  Energy  Transfer 

Energy  Flow  Rate 

1. 

impinging  water 

Ww  iWir  As 

2. 

evaporation 

We  iv.sur  As 

3. 

water  flow  out  of  control  volume 

Wr  *w,  sur 

4. 

ice  accumulation 

Wj  iijSur  As 

5. 

convection 

Qc  As 

6. 

conduction  through  skin 

Qk  As 

Using  the  convention  that  energy  flow  into  the  control  volume  is  positive,  the  terms  can 
be  summed  to  yield  the  energy  equation. 

Ww  iWjT  As  =  We  iViSU,  As  +  wr  iW|Sur  +  Wj  ii  sur  As  +  qc  As  +  qK  As  (B-4) 

The  evaluation  of  the  terms  of  the  energy  equation  has  been  done  by  various  authors, 
most  notably  by  Sogin  (Ref.  B-l),  Lowzowski  et.  al.  (Ref.  B-2),  and  Cansdale  and  Gent 
(Ref.  B-3).  The  following  sections  will  evaluate  each  of  the  terms,  highlighting  differences 
from  previous  models.  The  references  should  be  consulted  for  further  details. 

IMPINGING  WATER 

Since  the  droplets  are  essentially  brought  to  rest  when  they  strike  an  object,  it  is 
appropriate  to  use  the  stagnation  enthalpy  defined  as 

ul, 

iw.T  —  Cpw,s  (Ts  —  32)  +  (B-5) 

The  arbitrary  reference  for  zero  enthalpy  used  in  this  study  is  water  at  32°F,  Substituting  Eq. 
(B-5),  the  energy  flow  rate  of  the  impinging  water  becomes 

Ww  iw,r  As  =  Ww  I  ^Pw,s  —  32)  +  — — 1  As  (B-6) 

L  2gcJJ 

EVAPORATION 

The  rate  of  energy  transfer  from  the  surface  because  of  evaporation  is  given  by 
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we  i».sur  As  -  We  [cPw  sur  (Tsur  -  32)  +  Lv]  As  (B-7) 

where  We  is  the  evaporative  mass  transfer  rate  and  Lv  is  the  latent  heat  of  vaporization. 


The  mass  transfer  rate  is  analogous  to  the  convective  heat-transfer  rate  and  can  be 
written  as 


We  =  gAB 


(B-8) 


where  g  is  the  mass  transfer  coefficient  and  AB  is  the  evaporative  driving  potential.  The  mass 
transfer  coefficient,  g,  can  be  evaluated  using  the  analogy  to  heat  transfer  given  by  the 
equation 


(B-9) 


where  Pr  is  the  Prandtl  number  and  Sc  is  the  Schmidt  number. 


The  mass  transfer  driving  potential  is  analogous  to  the  temperature  difference  in  the 
convective  heat-transfer  equation.  In  the  case  of  evaporation,  the  driving  potential  is  a  vapor 
concentration  difference  instead  of  a  temperature  difference.  The  equation  used  in  this 
study  is  similar  to  that  derived  by  Sogin  (Ref.  B-l)  and  given  as 

ah  -  pv.*/Ts  -  <*VTt)  Pv>s/Ps 

(1/0.622)  Pt/Tt  -  PVlS/Ts  (B_I°) 

This  term  accounts  for  compressibility  effects  as  does  the  term  derived  by  Cansdale  in  Ref 
B-3. 


WATER  FLOW  OUT  OF  THE  CONTROL  VOLUME 

The  water  flowing  out  of  the  control  volume  will  be  at  the  surface  temperature,  allowing 
the  enthalpy  to  be  expressed  as 

Vsur  =  Cp^suf  (Tsur  -  32)  (B-l  1) 

Using  Eq.  (B-3),  the  runback  water  energy  flow  rate  can  be  expressed  as 

wr  'w.sur  =  [(1  -  n)  WwAs  -  WeAs]cPw  sur(Tsur  -  32)  (B-l 2) 
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ENERGY  OF  THE  ICE  ACCUMULATION 

From  the  definition  of  the  freezing  fraction,  Eq.  (B-2),  the  freezing  rate  is 


WjAs  =  n  Ww  As 


(B-13) 


The  enthalpy  of  ice  referenced  to  water  is  32°  F  is 

■i.sur  ~  cpiiSMr  l*^siir  “  32)  —  Lp  (B-14) 

Combining  Eq.  (B-13)  into  Eq.  (B-14),  the  energy  slorage  rale  in  the  control  volume  can  be 
expressed  as 


'i.sur  As  —  nWw  [^pj%sl|r(TM„  32)  —  Lj]  As  (It- 15) 

NET  CONVECTIVE  HEAT  FLUX 


Convective  heat  flow  rale  is  normally  defined  by  the  following  equation: 


q,.As=htATAs  (B-16) 

In  an  icing  analysis,  AT  is  the  difference  between  the  surface  temperature  and 
wall  temperature  and  can  be  expressed  by  following  equation: 

AT  =  T  -  T 

“ 1  1  xur  4  aw 

The  adiabatic  wall  temperature,  Taw,  is  given  by 

Taw  =  Tj  ^1  ■+  re  — - mA  (B-18) 


the  adiabatic 

(B-17) 


where  Tl  and  Ml  are  the  temperature  and  Mach  number  at  the  edge  of  the  boundary  layer, 
respectively,  and  rc  is  the  recovery  factor.  At  the  stagnation  point,  Ty  is  equal  to  the  total 
temperature,  Ty,  and  My  —  0.0.  The  total  temperature  can  be  expressed  in  terms  of  the 
free-stream  velocity  as 


T,  =  Tjj  + 


Uf 


Zfic-ltp., 


(B-19) 
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Substituting  Eqs.  (B-17),  (B-18),  and  (B-19)  into  Eq.  (B-16),  the  expression  for  the 
convection  heat  flow  rate  becomes 


Qc  As  =  hc 


-  Ts 


-  As 

2gc-lCpa  - 


(B-20) 


The  heat-transfer  coefficient  for  the  stagnation  line  of  a  cylinder  and  the  leading  edge  of 
airfoils  is  given  in  Ref.  B-4  as 


1.14  (Re)°-S  (Pr)0  4  k 
hc  - - 

c 


(B-21) 


where  Re  is  the  free-stream  Reynolds  number  and  c  is  the  cylinder  diameter  or,  in  the  case  of 
an  airfoil,  the  leading-edge  diameter. 


A  more  general  derivation  of  this  term  is  given  in  Ref.  B-2 
CONDUCTION  INTO  THE  SKIN 


When  the  cloud  is  first  encountered,  a  temperature  difference  will  exist  between  the 
wetted  surface  and  the  inner  structure  of  the  skin  which  was  at  an  equilibrium  temperature 
prior  to  entering  the  cloud.  The  evaluation  of  the  resulting  conductive  heat  flow  rate  is 
dependent  on  knowing  the  thermal  conductivity  and  detailed  geometry  of  the  skin.  The 
calculation  of  this  heat  flow  rate  through  what  could  be  a  complicated  composite  structure  is 
avoided  in  scaling  studies  by  requiring  that  the  objects  have  the  same  internal  strucure  and 
are  made  of  identical  materials.  It  is  then  assumed  that  the  conductive  heat  flow  rate 
between  the  model  and  full-scale  will  be  identical,  i.e. 


Qk  (As)model  Qk  (As)fu[]  -scale 

To  simplify  the  solution  of  the  energy  equation,  qKAs  is  dropped  from  the  equation.  This 
will  not  change  the  values  of  any  of  the  calculated  similitude  conditions  but  may  effect  the 
absolute  value  of  the  calculated  surface  temperature  and  freezing  fraction. 

The  energy  terms  can  now  be  summed  to  form  the  complete  energy  balance  equation  for 
a  stagnation-line  control  volume. 
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Ww  [c^CTs  -32)  +  2^j  J  "  We  j^CPw,su/Tsur  32)  +  LVJ 

+  [0  -  n)  Ww  -  WJ  cPw  itur  (Tsur  -  32)  +  n  Ww  [cPisur(Tsur  -  32)  -  Lf]  (B-22) 

+  he  [t,„  -  T,  -  ^-] 

Note  that  the  control  volume  surface  length,  As,  has  been  cancelled  uniformly  from  each 
term  of  the  energy  equation. 
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Figure  B-l.  Identification  of  a  stagnation-line  control  volume. 
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1. 

Impinging  Water 

Ww 

A  S 

2. 

Evaporation, 

W 

e 

AS 

3. 

Water  Flow  Out  of  CV 
(Runback) 

w 

r 

4. 

Ice  Accumulation 

Wi 

As 

Figure  B-2.  Mass  balance  for  a  stagnation-line  control  volume. 


1.  Impinging  Water  lw  T  As 

2.  Evaporation  le  ivaur  is 

3.  Water  Flow  out  of  CV  wr  lw>  auP 

4.  Ice  Accumulation  Wt  14  gur  As 

5.  Convection  q  As 

c 

6.  Conduction  As 


Figure  B-3.  Energy  balance  for  a  stagnation-line  control  volume. 
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NOMENCLATURE 

Ac  Accumulation  parameter 

AB  Mass  transfer  driving  potential 

b  Relative  heat  Factor 

Cd  Drag  coefficient 

c  Characteristic  length 

Cp  Specific  heat,  Btu/lbm 

D  Drag  force,  lbf 

d  Droplet  diameter,  ft  unless  specified 

dm  Mass  median  droplet  diameter,  ft  unless  specified 

Em  Total  collection  efficiency 

F  Equivalent  force,  lbf 

g  Mass  transfer  coefficient,  lbm-°F/Btu 

gc  Constant,  32.174  Ibm-ft/Ibf-sec2 

H  Projected  height,  ft 

h  Convective  heat-transfer  coefficient,  Btu/ft2/hr/°F 

i  Enthalpy,  Btu/lbm 

J  Mechanical  equivalent  of  heat,  778  ft-lbf/Btu 

K  Inertia  parameter 

Kq  Modified  inertia  parameter 

K  Particle  scaling  parameter 
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k  Thermal  conductivity,  Btu/ft/hr/°F 

Lf  Heat  of  fusion,  Btu/lbm 

Lv  Heat  of  vaporization,  Btu/lbm 

LWC  Liquid-water  content,  g/m3 

M  Mach  number;  Moment,  ft-lbf 

m  Mass,  Ibm 

m'  ‘  Mass  per  unit  area,  lbm/ft2 

n  Freezing  fraction 

P  Pressure,  psia 

Pr  Prandtl  number 

q  Dynamic  pressure,  psia 

qc  Convective  heat  flux,  Btu/ft2/hr 

qK  Conductive  heat  flux,  Btu/ftVhr 

R  Droplet  Reynolds  number  based  on  |  Ua  -  U<j  | 

Re  Free-stream  Reynolds  number 

Ru  Droplet  Reynolds  number  based  on  Ua 

St  Reynolds  number  integration  variable 

r  Droplet  radius,  ft 

s  Surface  length,  ft 

As  Surface  length  increment,  ft 
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Sc  Schmidt  number 

T  Temperature,  °F 

t  Time,  sec 

U.  Velocity,  fps 

W  Mass  flux,  lbm/ft2/sec 

w  Mass  flow  rate,  lbm/sec 

x  Length  integration  variable,  ft 

yo  Vertical  distance  between  tangent  trajectory  release  points,  ft 

0  Local  collection  efficiency 

y  Slope  of  log  CD  versus  log  R  curve;  ratio  of  the  specific  heats  of  air  =  1 .4 

0  Air  energy  transfer  potential,  °F 

X/Xj  Droplet  range  parameter 

P  Viscosity,  lbf-sec/ft2 

q  Density,  lbm/ft3 

ff  Shear  stress,  lbf/ft2 

r  Icing  time,  sec  unless  specified 

0  Droplet  energy  transfer  potential,  °F 

Subscripts 
a  Air 

aw  Adiabatic  wall 
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c  Critical 

d  Droplet 

e  Evaporation,  evaporative 

i  Ice 

1  Lower  surface 

L  Condition  at  the  edge  of  the  boundary  layer 

r  Runback  water 

s  Static  condition 

sur  Surface  condition 

T  Total  condition 

u  Upper  surface 

v  Vapor 

w  Liquid  water 

oo  Free -stream  condition 

Superscript 

'  Nondimensionalized  parameter 

—  Vector  quantity 
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